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Введение. Исследованы поперечные колебания биморфа, состоящего из двух пьезомагнитоэлектрических слоев 
и находящегося в переменном магнитном поле. Пьезомагнитоэлектрические слои представляют собой 
многослойный композит с чередующимися пьезоэлектрическими и пьезомагнитными слоями. Механические и 
физические свойства такого композита задаются известными эффективными константами. 

Материалы и методы. Прикладная теория колебаний многослойной пластины учитывает нелинейное 
распределение электрического и магнитного потенциала в пьезоактивных слоях в продольном и поперечном 
направлениях. На основе указанной теории исследованы напряженно-деформированное состояние, зависимости 
прогиба, электрического и магнитного потенциалов от объемного соотношения состава шарнирно опертого 
биморфа. Электрический потенциал принят равным нулю на всех электродах, магнитный равен нулю на 
внутренней границе и неизвестен на внешних. Поэтому распределение электрического и магнитного 
потенциалов в середине слоя — неизвестные функции. В случае магнитного потенциала распределение на 
внешней границе также неизвестно. В задаче были приняты гипотезы Кирхгофа для механических 
характеристик. Использование вариационного принципа и квадратичной зависимости электрического и 
магнитного потенциалов по толщине пьезоактивных слоев позволило получить систему дифференциальных 
уравнений и граничных условий. 

Результаты исследования. При изменении объемного соотношения состава пьезоактивных материалов 
биморфа электрический потенциал в середине слоя нелинейно изменяется. Магнитный потенциал в середине 
слоя и на внешней границе практически линейно увеличивается с увеличением объемного процента ВаТО:. 
Определена зависимость прогиба в середине слоя. 

Обсуждение и заключения. Построена прикладная теория расчета поперечных колебаний биморфа с двумя 
пьезомагнитоэлектрическими слоями. Исследована зависимость характеристик напряженно-деформированного 
состояния, электрического и магнитного полей от объемных долей пьезомагнитного и пьезоэлектрического 


материалов. 


Ключевые слова: пьезоэлектрики, пьезомагнетики, композит, биморф, магнитоэлектроупругость, изгибные 


колебания. 
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Пьезоэлектрические материалы напрямую преобразуют электрическую энергию в механическую и 
обратно. Это свойство позволяет широко использовать их в науке и технике. Данные материалы задействуют в 
ультразвуковых излучателях упругих и акустических волн, приемниках таких волн, устройствах подавления 
колебаний элементов машин и конструкций и т. п. В последнее время бурно развивается еще одна область 
применения пьезоэлектриков — устройства сбора и накопления энергии. В этом случае пьезоэлектрические 
материалы входят в состав пьезоэлектрических генераторов энергии (ПЭГ). ПЭГ помещаются на элементах 
машин или конструкций, которые интенсивно колеблются, находятся в зоне прохождения упругих волн или 
подвержены действию переменного давления. Основные типы этих устройств имеют биморфную или стековую 
многослойную структуру и испытывают изгибные или продольные деформации соответственно. На основе 
ПЭГ создаются маломощные источники электрического тока. В их числе — автономные источники питания 
(например, для устройств мониторинга повреждений в труднодоступных местах сооружений трубопроводов 
ит. п.). Обзор таких устройств есть в [1-2]. Один из способов конструировать эффективные ПЭГ — 
использование пьезоактивных композитов различного типа связности и неоднородных материалов на основе 
пьезокерамики, в том числе пористой. 

Фиксировать или задействовать энергию переменного магнитного поля могут ПЭГ, в конструкции 
которых есть дополнительные электромагнитные элементы или постоянные магниты. Один из путей решения 
данной задачи — использование пьезомагнитных материалов в сочетании с пьезоэлектрическими. В этом 
случае переменное магнитное поле приводит к деформации пьезомагнетика и связанного пьезоэлектрика, в 
результате последний генерирует электрическую энергию. Существует класс материалов, обладающих 
пьезомагнитными свойствами. Пьезомагнетизм — это явление, наблюдаемое в некоторых антиферромагнитных 
и ферромагнитных кристаллах. Он характеризуется линейной связью между магнитной поляризацией системы 
и механической деформацией. В пьезомагнитном материале можно вызвать спонтанный магнитный момент, 
приложив механическое напряжение, или деформацию, приложив магнитное поле. В исследованиях 
пьезомагнитных материалов очень часто рассматривается СоРе›О. [3-5]. В [6-8] исследуется композит на 
основе СоГе>О; и ВаТ1Оз, обладающий пьезоэлектрическими и пьезомагнитными свойствами одновременно. 

Решения задач электроупругости и магнитоупругости приведены в [9—11]. В [12] развиты прикладные 
теории колебаний многослойных пьезоэлектрических пластин с учетом специфики распределения 
электрического потенциала по толщине конструкции. 

Задачи об установившихся колебаниях электромагнитоупругого слоя и полупространства под 
действием гармонических нагрузок представлены в [13, 14]. Учтены предварительные напряжения, а также 
различные электрические и магнитные условия на границах. Исследовано влияние указанных факторов на 
дисперсионные свойства. 

Ранее [15, 16] была разработана прикладная теория, которая учитывает неоднородное распределение 
электрического потенциала в продольном направлении и квадратичную зависимость по толщине. В этих же 
работах исследовано напряженно-деформированное и электрическое состояние шарнирно опертого и 
консольно закрепленного биморфа. В обоих случаях прикладная теория показала хорошую сходимость с 
результатами конечноэлементного моделирования. Авторы также получили прикладную теорию колебаний 
биморфа [17], состоящего из электроупругого и магнитоупругого слоя. Этот подход хорошо согласуется с 
результатами конечноэлементного анализа. 

В данной работе рассматриваются колебания устройства в рамках плоской деформации. На основе 
вариационного принципа построена прикладная теория изгибных колебаний — двухслойного 
пьезомагнитоэлектрического биморфа. Для установившихся колебаний получены граничные условия и система 
дифференциальных уравнений для четырех неизвестных функций (прогиб, электрический потенциал в 
середине слоя, магнитный потенциал в середине слоя и на внешней границе), зависящих от длины биморфа. 
Исследовано влияние разных процентных объемных соотношений состава биморфа на прогиб, электрический и 
магнитный потенциалы в определенных положениях. Результаты исследования позволяют выбрать состав 
композиционного пьезомагнитоэлектрического материала для достижения наиболее эффективной работы 
устройства. 

Материалы и методы. Рассматривается пластина, состоящая из двух одинаковых 
пьезомагнитоэлектрических слоев. Она совершает установившиеся поперечные колебания в рамках плоской 
деформации. Каждый слой представляет собой композит 2—2 связности, состоящий из чередующихся 
пьезоэлектрических и пьезомагнитных слоев (рис. 1). 
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иашиии  |[ьезоэлектрический слой м — Пьезомагнитный слой 
Рис. 1. Структура композита 2—2 


В [8] найдены эффективные свойства такого композита. Большие поверхности слоев электродированы, 
а сами слои поляризованы по толщине. Биморф шарнирно закреплен по краям, все поверхности свободны от 
механических напряжений. На верхнюю и нижнюю границы пластины воздействует магнитный поток Ву, в то 
время как на границе между слоями магнитный потенциал считается равным нулю. Электрический потенциал 
равен нулю на всех электродах. Боковые поверхности считаются изолированными от магнитных и 
электрических полей. 
Уравнения для описания колебаний композита с эффективными свойствами, связностью механических, 
электрических и магнитных полей, имеют вид [18]: 
У.о+р/ =рй, У.-Р=о., У:В=0, 
о=с:=-е'.Е-Й’.Н, 
Р=е:=+к.Е+а-Н, 
В=Й:=+а' .Е+и-Н, 


в=5(уи+ (и) Е=-Уф, Н=-УЕ. (1) 


Здесь би & — тензоры механических напряжений и деформации, ОД и Е — векторы электрической индукции 
и напряженности электрического поля, В и Н — векторы магнитной индукции и напряженности магнитного 


поля, р — плотность материала, с — тензор упругих модулей, е — тензор пьезоэлектрических модулей, р — 


тензор пьезомагнитных модулей, к —[ тензор диэлектрических проницаемостей, [0 —[_ тензор 
магнитоэлектрических модулей, и — тензор магнитных проницаемостей, Л —Щ_ вектор плотности массовых 
Сил, бо — объемная плотность электрических зарядов, и — вектор перемещений, фи |2 == электрический и 


магнитный потенциалы. 

Граничные условия определяются для механического, электрического и магнитного Поля 
соответственно. 

Для первого случая отметим отсутствие механических напряжений на границе биморфа: 


уп]. =0, &1=13. 


Биморф шарнирно закреплен на концах (рис. 2): 
= 0, 1=1,3. 





и о = и 
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Рис. 2. Геометрия и граничные условия биморфа с композитными пьезомагнитоэлектрическими слоями 


Далее сформулируем электрические граничные условия. Электрический потенциал на внутреннем и на 
внешнем электроде соответственно: 


Е = и, Фи = , у 


Укажем магнитные граничные условия. Магнитный потенциал на внутренней границе: 
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6 йе М.. 
Магнитный поток Во воздействует на верхнюю и нижнюю границы пластины: 
Ж=ЕН = Ву у 


Воспользуемся вариационным уравнением для установившихся колебаний [10]. Оно обобщает принцип 
Гамильтона в теории электроупругости с учетом магнитных составляющих. Для случая плоской деформации 
при отсутствии поверхностных нагрузок и при наличии магнитного потока: 


[584$ -роу [[изби,а$ + [| ( р,ди, +в,6ф+ В,6Е)4$ =0, (2) 
К 5 5 


где ЭН = с,6=, -О5Е, — ВН, . 
Для построения прикладной теории колебаний примем гипотезы Кирхгофа. В соответствии с ними 
распределение перемещений по толщине имеет вид: 
ми ьх,) = хуи, ь и,(х,х,) = их) (3) 


В частности, для механического поля принята гипотеза единой нормали. Далее рассматривается задача, 
В которой значение электрического потенциала на электродах может быть равно нулю, поэтому его 
распределение не описывается линейной функцией. С учетом возможной неоднородности по длине элемента, 
связанной с влиянием граничных условий на концах биморфа, его распределение по толщине принимается 


Ф(х,5)=И,(х) 22 | ксу(- | и<)-> 228) (4) 


Здесь х, =х,-Н/2. Функции Г,Г и Г, отвечают за значение электрического потенциала на внутреннем 


квадратичным: 





1 
электроде, в середине слоя и на внешнем электроде соответственно. Чтобы удовлетворить условиям задачи, 
примем эти функции в следующем виде (см. рис. 2): 
И.) =И, =соп5ь И (х)=Ф(х,), Г, (2%) ЕТ, = соп&1. 
Здесь функция Ф(х,) является неизвестной. 


Представим квадратичное распределение магнитного потенциала по толщине каждого слоя. 
Распределение по длине неоднородно, на внутренней границе слоев его значение принимается равным нулю: 


оу ми ме су +М, (=) Е =) (5) 


Здесь Хх, =х-Н/2. Функции М, М, и М, отвечают за значение магнитного потенциала на внутренней 
границе, в середине слоя и на внешней границе соответственно и принимаются в следующем виде (рис. 2): 
Мь(х) =М, =сопя, М, (х,)==.(х), М‚(х)==.(х). 
Здесь функции =,(х,) И Я, (х,) являются неизвестными. 
Соотношения (3)—(5) подставим в уравнение (2) и проинтегрируем его по толщине биморфа, а затем 


приравняем к нулю коэффициенты при независимых вариациях би, БФ, б=, и бЕ, . Тем самым мы получаем 


систему из четырех дифференциальных уравнений (6) от четырех неизвестных функций, зависящих от Хх, 
(далее опустим нижний индекс), и пять граничных условий (7). 
1 6Е. 16Е. 32Е. 32а 164 : 
633 + 633 и 63 Ф(л) ЕЕ =,(х)+ 64. =, (х) 16. Н — Ф(^) 
ЗН ЗН ^ ЗН ЗН - ЗН 15 4 
_ 16а, Н 4” 2а,Н а? 46, Н а? 16а 
ЕВЕ 
о о к” 
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т И у. = Ф(5) = = те 5 вый е а 
ре г =.) та те а ый Е ко м 
Ру, м + т Ф)+ Е ©) М =.) - е те 

Вы, ми 5,69 + Зы 4 О-о, -2вм, 0 
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4г,Н — 4йН а? 5й,Н 4? 
—— Ф(х) + ы На ) ы 523 (х) 2р-Н 
3 3 3 м (6) 
НН? а? Е 
р #4” с, Ч сих =0. 
м 3 @ 
16Е.Н 1 Н 2а,Н 
т о Е Е 50 
15 К 15 @ 15 
160 Н а 164, а 2н:Н а — 
Ф( хх) + = =. (х)=0, 
5 15 
20 Н а 2ниН а — Ан а — 
Ф(х) + =, (х — =. (х)=0, 
а в 
НН, 5е,Н., 4е.Н 4йН _ 5й.Н — ЭН", 9" НЫ 
| 7 + Ф(х я. (х =. (х)-+ —=и(х)+ М, =0, 
в и (+ : 2(%) - 3(%) - а, о 
46, 4 4#,Н а 5й.Н 4 — НЕ. а? 2ро?Н? 4 
Ф я, (х)+ =. (х их) -+ их) = 0. 7 
те (х) - Е. ты 5 аа ие (7) 
Здесь были введены следующие — обозначения: а=а=а Ла. &=е ее 0. 
Й, =й, -с,й, /с,, @, =—а,-е,й,/с,, Е, =-е, —е. /с.. Они возникли после удовлетворения условию 


„ =0 и исключения 5... 


Результаты исследования. Сопоставлены результаты расчета биморфа по предложенной теории с 
конечноэлементным расчетом в низкочастотной области для соотношения объемных долей 
пьезоэлектрического и пьезомагнитного компонентов 80 % ВаТ.Оз и 20 % СоЕе>О4л. Сравнение показало, что 
погрешность в нахождении характеристик механического и магнитного полей менее | %. При определении 
электрического поля в средней части пластины разница составила порядка 5 %. Описывая ситуацию в 
окрестности точек опоры, следует отметить, что размер окрестности по продольной координате 
приблизительно равен толщине биморфа. Здесь при определении электрического поля зафиксирована 
разница 20 %. 

Первым шагом исследования колебания двухслойного пьезомагнитоэлектрического биморфа при 
изменении объемного соотношения Ва; и СоЁе›О4 в составе композита является определение его 
эффективных свойств. В таблицах | и 2 представлены эти свойства, найденные по результатам работы [8]. 












































Таблица 1 
Материальные константы (упругие модули, диэлектрические и магнитные проницаемости) 
для различной объемной доли ВаГОз 
Объемная ен Диэлектрические Магнитные 
доля проницаемости проницаемости 
Ват: ГПа 10’ Ф/м 10° Нс/Кл” 
(%) Си С С1з Сзз Сад К Кзз И Изз 
0 286,0 173,0 170,0 269,5 45,30 0,080 0,093 5,900 1,570 
10 270,9 160,4 154,9 260,0 45,07 1,469 0,073 5,315 0,632 
20 256,6 148,5 142,6 250,2 44,84 2,815 0,098 4,730 0,396 
30 242,8 137,2 131,3 240,8 44,61 4,063 0,122 4,145 0,285 
40 229,9 126,8 120,9 231,9 44,38 5,287 0,147 3,560 0,223 
50 217,6 116,9 111,0 224,0 44,15 6,413 0,171 2,975 0,186 
60 206,7 108,1 102,1 215,6 43,92 7,490 0,220 2,390 0,155 
70 195,9 99,7 93,8 208,2 43,69 8,517 0,294 1,805 0,136 
80 186,0 92,3 85,9 201,3 43,46 9,448 0,441 1,220 0,120 
90 176,6 85,4 78,9 193,9 43,23 10,353 0,857 0,635 0,110 
100 166,0 77,0 78,0 162,0 43,00 11,200 12,600 0,050 0,100 
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Таблица 2 
Материальные константы (пьезоэлектрические, пьезомагнитные и магнитоэлектрические модули) для 
различной объемной доли ВаТ!Оз 












































Объемная Пьезоэлектрические модули Пьезомагнитные модули Магнитоэлектрические модули 
доля Кл/м” Н/А м 10° Не/ВКл | 10Т'Не/ВКл 
ВаТЮз 
(%) @1 @зз 15 п. п. 5 и б..; 
0 0 0 0 580,3 —699,7 550 0 0 
10 —0,006 0,029 1,16 223,6 —244,1 495 —1,33 1,97 
20 —0,013 0,059 2.32. 130,0 —132,3 440 —2,35 2,36 
30 —0,019 0,088 3,48 86,7 —79,8 385 —3,07 2,48 
40 —0,025 0,132 4,64 61,6 —52,5 330 —3,48 2,50 
50 —0,031 0,176 5,80 43,3 —34,2 275 —3,62 2,47 
60 —0,038 0,220 6,96 29,7 —22,8 220 —3,45 2,43 
70 —0,040 0,352 8,12 20,5 —13,7 165 —3,00 2,36 
80 —0,060 0,571 9,28 13,7 —9,1 110 —2,27 2,29 
90 —0,263 1,187 10,44 4,6 —4,6 55 —1,28 2,16 
100 —4,400 18,600 11,60 0 0 0 0 0 



































Колебания биморфа возбуждались магнитным потоком, приложенным к верхней и нижней 
граням (рис. 2), который изменялся по гармоническому закону с амплитудой В, =5х10® Вб и частотой 10 кГц. 

На рис. 3 представлен прогиб в середине слоя в зависимости от объемной доли Ваз. Из графика 
видно, что прогиб в положении, имеющем координаты х =[/2, х, =Н/2, равен нулю, если биморф состоит 
только из пьезоэлектрика ВаТ!Оз. Прогиб биморфа достигает наибольшего значения, если в его составе только 
пьезомагнетик СоЁе›О.. Прогиб практически линейно зависит от объемного соотношения составляющих 
пьезоактивных материалов. 

Прогиб им (),м х 107" 
00 


0 20 40 60 80 100 
Объемная доля ВаТ!Оз, % 


Рис. 3. Прогиб и(х,) в середине слоя для различной объемной доли Ваз 


На основе данных рис. 4 можно сделать заключение, что электрический потенциал в середине слоя 
нелинейно изменяется при изменении объемного соотношения состава пьезоактивных материалов биморфа. 
Если биморф состоит только из ВаГЮз либо СоРе>О%, то электрический потенциал в точке (2/2,Н/2) равен 


нулю и достигает наибольшего значения при 35 % ВаТЕО: в составе биморфа. 
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Рис. 4. Электрический потенциал Ч(х,) для различной объемной доли Ва1Оз 


Анализ рис. 5 и 6 позволяет сделать заключение, что магнитный потенциал в середине слоя =,(Ё/2) и 


на внешней границе =,(Ё/2) практически линейно увеличивается с ростом объема ВаТ1О: в составе биморфа. 
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Рис. 5. Магнитный потенциал =,([/2) для различной объемной доли Ваз 
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Рис. 6. Магнитный потенциал #.(Ё/2) для различной объемной доли Ваз 


Обсуждение и заключения. Предложена прикладная теория для расчета поперечных колебаний 
биморфа из двух слоев композита на основе СоРе›О. и ВаПГОз, обладающего одновременно 
пьезоэлектрическими и пьезомагнитными свойствами, в переменном магнитном поле. Такая конструкция И 
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может служить моделью пьезоэлектрического генератора устройства сбора и накопления энергии при действии 
внешнего магнитного поля. В низкочастотной области (ниже собственной частоты первой изгибной моды) 
проведены расчеты напряженно-деформированного состояния биморфа, распределения электрического и 
магнитного полей. Исследована зависимость прогиба, электрического и магнитного потенциалов от объемного 
соотношения состава биморфа. В дальнейшей работе предполагается определить выходной потенциал и 
мощность электрического тока, возбуждаемого переменным магнитным полем. Целью этих изысканий будет 
сбор электрической энергии. 
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Введение. В статье изучается осесимметричная задача теории упругости для радиально-неоднородной 
трансверсально-изотропной незамкнутой сферы, не содержащей ни один из полюсов 0 и п. Считается, что 
модули упругости являются линейными функциями от радиуса сферы. Предполагается, что боковая 
поверхность сферы закреплена, а на конических сечениях заданы произвольные напряжения, оставляющие 
сферу в равновесии. 

Целью данной работы является асимптотический анализ задачи теории упругости для радиально-неоднородной 
трансверсально-изотропной сферы малой толщины и исследование на базе этого анализа трехмерного 
напряженно-деформированного состояния. 

Материалы и методы. На основе уравнений теории упругости методом однородных решений и 
асимптотического анализа исследуется трехмерное напряженно-деформированное состояние радиально- 
неоднородной сферы. 

Результаты исследования. После выполнения однородных граничных условий, заданных на боковых 
поверхностях сферы, получено характеристическое уравнение и произведена классификация его корней 
относительно малого параметра, характеризующего толщину сферы. Построены соответствующие 
асимптотические решения, зависящие от корней характеристического уравнения. Показано, что решения, 
соответствующие счетному множеству корней, имеют характер пограничного слоя, локализованного в 
конических срезах. Разветвление корней порождает новые решения, которые характерны только для 
трансверсально-изотропной радиально-неоднородной сферы. Появляется слабозатухающее погранслойное 
решение, которое может проникать глубоко вдали от конических сечений и изменять картину напряженно- 
деформированного состояния. 

Обсуждение и заключения. На основе построенных решений можно определить области применимости 
существующих прикладных теорий и предложить новую более уточненную прикладную теорию для радиально- 
неоднородной трансверсально-изотропной сферической оболочки. 


Ключевые слова: уравнения равновесия, уравнения МЛежандра, радиально-неоднородная сфера, 
характеристические уравнение, погранслойные решения, вариационный принцип, прикладная теория, метод 


редукции. 
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Введение. Одним из свойств материалов, влияющих на напряженно-деформированное состояние 
упругих тел, является их неоднородность. Изучение напряженно-деформированного состояния неоднородных 
тел на основе трехмерных уравнений теории упругости связано со значительными математическими 
трудностями. 

Исследованию трехмерных задач теории упругости для сферы посвящен ряд исследований. 

В работе [1] на основе уравнений теории упругости для сферы получено общее решение, 
удовлетворяющее граничным условиям на контуре в смысле Сен-Венана, проведен анализ напряженно- 
деформированного состояния сферы. В [2] на основе уравнений теории упругости для толстой изотропной 
сферы построены однородные решения, зависящие от корней трансцендентного уравнения. В [3] на основе 
решения трехмерных задач теории упругости для сферы малой толщины изучена точность существующих 
прикладных теорий и дан метод построения уточненных прикладных теорий. В [4] изложена трехмерная 
асимптотическая теория трансверсально-изотропной сферической оболочки малой толщины. В [5] приведен 
анализ трехмерного напряженно-деформированного состояния трехслойной сферы с мягким заполнителем. 
В [6] методом однородных решений изучена задача кручения для радиально-неоднородной трансверсально- 
изотропной сферы малой толщины, когда упругие характеристики меняются линейным, квадратичным и 
обратно квадратичным законами по радиусу. В [7] изучена задача кручения для радиально-слоистой сферы с 
произвольным числом чередующихся жестких и мягких слоев. Показано существование слабозатухающих 
погранслойных решений и возможное нарушение принципа Сен-Венана в его классической формулировке. 
Построена прикладная теория кручения радиально-слоистой сферы, адекватно учитывающая возникающие 
особенности. В [8] с помощью метода конечных элементов и сплайн-коллокации исследована задача теории 
упругости для радиально-неоднородного полого шара. Проведено сравнение результатов, полученных 
методами конечных элементов и сплайн-коллокации. В [9] методом асимптотического интегрирования 
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уравнений теории упругости изучена осесимметричная задача теории упругости для радиально-неоднородной 
трансверсально-изотропной сферы малой толщины. Построены неоднородные и однородные решения. 
Установлен характер напряженно-деформированного состояния. В [10] методом однородных решений 
рассмотрена осесимметричная задача теории упругости для сферы малой толщины с переменными модулями 
упругости. Получены асимптотические формулы для перемещений и напряжений, позволяющие рассчитать 
трехмерное напряженно-деформированное состояние радиально-неоднородной сферы. 

Материалы и методы. Рассматривается деформация в рамках линейной теории упругости 
незамкнутой сферы, материал которой является неоднородным по радиальной координате и трансверсально- 
изотропным. Толщина полой сферы предполагается малой, по сравнению с радиусом и размером по дуговой 
координате. Рассматриваются граничные условия, позволяющие решать задачу в осесимметричной постановке. 
Предполагаем, что сфера не содержит ни один из полюсов 0 и л. В сферической системе координат область, 


занятую сферой, обозначим через Г = {" =["; |], 9е[@; 60,], ФЕ[0; 27}. 
Рассматривается линейная зависимость упругих свойств материала по радиусу: 
(0) ©), — 70) = 
А =ал Г, А, = то Ной ‚ 4 =а> г 2 Аз; = аз Г, Ан =ад Г (1) 


(0) (0) (0) (0) (0) 
где, аи ›а> ›@>> ‚Аз › Ча 


—Щ_ некоторые постоянные величины. 
Система уравнений равновесия при отсутствии массовых сил в сферической системе координат г, 0,9 


имеет вид [11]: 





до, ь, 100% - 20. -0%-бъ +080 





=0, 2 

ОГ г 00 х _ 

00», 100 , 30. +9) 80 (3) 
О’ г 090 т 


где, с 


т? 


Оо Об > Обр —Щ_ компоненты тензора напряжении, которые выражаются через компоненты вектора 


перемещений о, =6, (г,0), о, =0,(т,0) следующим образом [4]: 








о.=А, = + А- ешо +20. + 7 , (5) 
о„=А, — +- Е (бы 5} о. + А›0, се + А,, 7 | (6) 
в = 4 : © +А,)0.+А.0,с1е0+ А, о | (1) 

си у 


Подставляя (5)-(8) в (2)—(3), с учетом (1) получаем уравнения равновесия в перемещениях. 











и ди 
ыр 52 "_ +2 (25% 5%?) =?и, +(250 5% 5% 5% }х 
ди, ди 0°и 0’и, би 
х=? | Очи себ |+ =(Ь®+Ь® ]| - бсюб-+—® +=? РР 40 |= 0, 9 
50 и г) (% и др 8 вер ч | 50? | 20 с (9) 


























д? д ^и ‚ ди 
Бо ы Пе т + (35 + +) 2+ 
др р дро@ 00 
21 (0) - (0) би, (0) (0) 
ч=?Ь® | в срб+ т += ве =? (Боев? 0 + и, =0. (10) 
1 а В и ь 1 7, 5 ы 
Здесь р=-Ш — новый безразмерный переменный; # = т. — | — малый параметр, характеризирующий 
8 п и 
(0) 
з ь 1. о, о _ Ч | 
толщину сферы; х=\ию; ре[-Б Ш и =, и, =—, В, = С безразмерные величины; ©, — 
^ ^ 0 


некоторый параметр, имеющий размерность модуля упругости. 
Предполагаем, что боковая часть границы сферы закреплена, т.е. 
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=0, (11) 
в: (12) 
Считаем, что на торцах сферы (на конических срезах) заданы напряжения 
оо = Л.Р), бы = Л.Р). (13) 
Здесь /,(р), №». (р) (5 =1;2) — достаточно гладкие функции, удовлетворяющие условиям равновесия. 
Решения (9), (10) ищем в виде [3, 4]: 
и,(р,0) =а(р)т(@); и„(р, в) =4(р)т"6), (14) 
где функция т(0) удовлетворяет уравнению Лежандра: 
итеучсшо-и\б)+[ 2 =" @ = 0. (15) 
После подстановки (14) в (9), (10), (11), (12) с учетом (15) получаем: 


Бру зараоунай (вы зы 28)-[2-1Ы? о)- 











(16) 
= - С +50) ар) = (+ +5 250 )а(р)|=0 
д 44 12 44 22 23 12 > 
В атоуеайр- [2 рр ар) р в 
=? (ЗЫ +5 +5 а(р)+=(5% +50) а(р)=0, 
а(р)=0, при р =-1. (18) 
а4(р)=0, при р=+1. (19) 
Решение системы (16), (17) имеет вид: 
ара" [ ре“, + ре С, + рае", + ве С, |, со) 
4(р)=е“ [ле С +4е “*С, чье” С, +4ье “*С я ь (21) 
где С, (п= 1,4) — произвольные постоянные, 
ры? [2-1] + (8 9-48); 
ь = (9 +557), — (255) +52 +53 -№5); 
а, = [5 ЕО = (2% Вр = ); 5; — корни уравнения 
нрырые [2-1 (у +знры вов) вры - ИРЫ «азы? - 
Ра ФО оО 2 + (2 - 1 4 +25 +29 + кр . 62) 


+ в 1 +50 50 +4 -(2 ь ЕН НЫ 0, 


Система линейных алгебраических уравнений относительно С,С,,С.,.С., получается путем 


удовлетворения однородным граничным условиям (18), (19). Равенство нулю определителя этой системы 
является условием существования ненулевых решений и приводит к характеристическому уравнению 
относительно спектрального параметра 2: 


4(2;=) —= (245 — 2.4, (и Р> -БР 57° (=(5 +5) )) + 
+( рр, — р,4, (рн, — рый) 5? (&(5, —5,))=0. 


Уравнение (23) имеет счетное множество корней 2,. Общее решение задачи получается путем 


(23) 


суммирования по корням уравнения (23) 


и, =УМ,а, (рут, (9), 4) 


= 


Механика 


17 


Бйр://уезиК-Чоп$а.га 


- 
со 


А4уапсей Епатеетия Кезеагсй 2022. Т. 22, № 1. С. 14-23. 155М№ 2687—1653 





=>. м4, (рт) вы 
где р 
брус [реа ето, ре (ета ео} 
а 


—=( +255) =(25.— *), 


+ 4 ы — Р>4! )е 





@ = р>4! (а. ь)е^ 1, (ра Р>4! )е 
©, = Р.& (а, -Ь а“ -Ь (р4> —Р.й )е 


=(я-—25 =(5 +25>) 
› 


г) +4› (Рь — Р.® )е 


—=(25, +55) 


—4, (214> — Р2й де в + Р19› (а, —& )е”* р 
@, =1 (рь — Р›4, мы — Ч (рь = 2в + Рь (а, —& )е = : 


Множество корней уравнения (23) при = —> 0 состоит из счетных множеств корней 


6, 


©, = (24. — Р.4, }е 


=( 251 +5> ) 











д =—и+0(е). (26) 
Для д, имеем: 
1°. При В >0, № -Ь, >0: 
(5—5) (6% 505.5, )з#((я +5,)6) (5, +5) (6% +605,5, эт ((5, -5,)6)=0, (27) 
где 
де; Д-р -Ь, 
в = (25°) в. +00 -5РЫР); В, = ОФ, 
2°. При В >0, Р-Ь, <0: 
выра? — (5 (а? - 8*)-5®) (26а) м 
ва [ 2508? + (65° (а В?) - 5?) | 5п(298)=0, 
где 
я= +62 -Ь, =(а +16); 
ЕЛЬ -Ь, =+(а-18) 
3°. При В >0, №=Ь: 
(5 (256) +2(675° +2) 65=0, (29) 


где 5 = в 
4°. При В <0, №-Ь,>0: 
(5, —5, (5% =, }5й ((5, +5,)6)+ (5 +5, (5% +505, }5й ((5-—5,)56)=0, (30) 
где 
я -Ь; зы Ь,. 
59. При В <0, 6? -Ь, <0: 
выра? (5 (22 - *)- 5 }| зи (26а) + 
(31) 
а 25 в? +(5° (а? °)-5° Е (258)=0, 


где 
х=ы- УР, +а-врх = ее -Ь, а 
60. При В <0, 2 =Ь: 
[В =6%) 51 (256)+2(Ы В) 5? +9) 65 =0 — 2) 
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где $=/Ь| 


Уравнение (27)—(32,) имеет счетное множество решений. 


Приведем асимптотическое построение решений, соответствующих различным группам корней 
характеристического уравнения (23). Подставляя (26) в (24), (25) и, раскладывая полученные выражения по 


степеням & , имеем: 
1 


а) и, (р;0)= УЕ 5 |(ы +5) (5% +5055 5, 605 (6,5, )зт (6,5 р)- 


КЕ 


(5% +552) 5, соз(6,5, )т (6,56) | +0(2) т, (6), 


(33) 


и (р;0)= УЕ‘ 5 (5 +5257 ) (5% +5055 | оз (6,5 )соз(5,5р)-с0$(5,5, )соз (6,552) ] +0(=)} т, (8), (34) 


= 


где д, являются решениями уравнения 


[8 5.) (5 -605 5, (5 +5,)5)+(5. +5, ) (5 +605,5, за (5. —5,)5)=0 


6) и,(р:0) =>. Е (5 +5) (5 +952) зп ( 645, )с0з(дяр)- 


= 


- (5% +5052 5, чп (буя )соз (дур) |+О(=)}т, (0), 


(35) 


(36) 


и, (р;0)= УЕ (52 +05 (5 +5755 )[ (9,5, т (б,5.р) -11 (6,5, 511 (6,5 р)]+0(=}} т, (6), (37) 


К=1 


где бд являются решениями уравнения 


(5, 5.) (6% =), 511 (5, +5)6)-(5 +5, (№ +155 511 (5, —5)6)=0. 
г. 
а) и, (р;0) = УЕ 5 (о +) {[ Вэ (5, 60) 51 (б,ар)- и соз(6,Вр)сй(6,ар)|. 
| авы? соз(5, 6) 51 (6,4) + (5 (в? — В*)-5® (5, В)ей (ба. | = 
-[а зт (5, р) 5й(9,@р)-+ В соз (,Вр)сй(б,ар)] . 


[рарыт чт (5,В)ен(ба) (5% (в? - В°)- 5) сов.) 51 (ба) |+ 0(е)} т(@), 


и, (р;6) = У Е 54 {авы зв (6, вр)св(днар)- (в (а? 1?) 5} 


КЕ 


жсоз(5, 6) 5й (6, @р) | [2авы® с0з(5,В)5й (6,&) +(5> (г -?)-5% 5 (5, В) ей(6,а) | - 


-[2авы? 0$ (68, р)5й (5, ар) ( (“г _ В* ) -5о т (6, Вр)сй (бар) | х 


х[2авы эт (6,в) ей (на) (5 (о? В?) --5Ф со (6,51 (буд) | +0(=)} т, (6), 


где д, являются решениями уравнения 


(38) 


(39) 


(40) 
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в[25тее — (5% (в? - в)- 5%) [58 (26а)+а[ 25? 8" + (5 (=-в?)-5$) (298) =0. (41) 


6) и, (р;0) = УЕ 5 [5 +5) {| В5т(5, Во) сй(6,ор)- и соз( 6, Вр) 51( 6, ар) |. 


| 2аВЫ? соз(6,В)ей(6,а)+ (5? (2 — В?)-5® 5 (6, В)5# (ба) |-- 
—[азт(6, Вр) сй(б,ар) + В соз(5, Вр) эй (б,ар) | 2аВ ВБ” зт(5,В) 51 (6,4) 
-(5 (2 -1°)- 55 )ооз(5, в) св( бла) + 0(ё)|т, (6), 


и, (р:6)=У`Е®5# [[орьо мп(5, Вр) 5й (бар) (5° (а? — в*)- 5х 


КЕ 


хеоз(5,Вр)св(б,ар) | 2аВЫ? соз(5,В)сй(6,а)+(5? (2 — °)--5$ 5 (5,8) 5й (5,2, |- 
-[2авы? со (58, р)сй(9,ар) +(5> (г — В ) -5 )5т(6,6р)5й (бар) | С 


х[2ары? зт(5,В)5й(5,&) =|8 (= _ 6? уе сз (д, В) ей (д.4 | +0(2)} т, (8), 


где д„, являются решениями уравнения 
в] 25 тее — (5? (в? — В°)-55) | #(26а) а 25% в? + (5? (2 -В°)--5$) (298) =0. 
—. 


а) и, (р;0)= УЕ (95? +0) (с03(6,5)с03(6,5р) + рз(6,5) зп(6,5р)) + 


К=1 


(0) ‚2 (0) 
(6 5 5. ) 


50 


т (5,5)с0$ (6,52) + 04) (8), 


2 
(о +59) 


и, (2;0)= > р 
х(р5й1 (5,5) с05(6,5р)- с0$(6,5)з(6,5р))+0(=)} т, (6), 
где бу являются решениями уравнения 
(67 5% (256) +2(6057 +5) 65=0 


В© 52 -5%) 
Е® ( и 44 
6) и, -> и - 


с0$(6,5)т(6,5р) - (55? +50) х 


х(51т (5,5) т (9,52) + рсоз(5,5) соз(5,5р)) + 0(=)} т, (9) 


(7 +) 
(6% +82)65 


и, (р;0) => Е | 


&= 


(42) 


(43) 


(44) 


(45) 


(46) 


(47) 


(48) 
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х| рсоз(6,5)зт (5,5р)-зт(9,5)с05(6,5р) |+ 0(=)} т, (в), (49) 


где д, являются решениями уравнения 
(7-5) т (256)-2(5757 +50) 65 =0 (50) 
4°. В случае В <0, Ь? -Ь, >0 асимитотические формулы для перемещений получаются из (33)-(38) 
заменой 51,5, на 15, ,15.. 
59. В случае Ь < 0, В? -Ь, <0 асимптотические формулы для перемещений получаются (39)(44) 
заменой 51,5, на 15,,15.. 
6°. В случае Ь < 0, 5? =Ь, все асимптотические формулы для перемещений получаются из (45)-(50) 


заменой 5на 15 


Для корней (26) главный член асимптотического решения уравнения (15) при =->0 принимает 
вид [9, 10]: 








- - яч ='\/-52 (0-6, ] (1+0(=)); в окрестности В =8, 
эт 4-92 
т, (в) = ь (© 
ке #288 (9-6, )| (1+0(=)); в окрестности @ =6,. 
тб + 5, . 
Перемещения представим в виде: 
и„(р,0) => Е,а, (рут, (6), (52) 
ри 
и (р, в) = > Е, а, (рт, (в). (53) 
К=1 
Напряжения о и о представим в виде: 
ош = УЕ, (о (рут, (в) +0 (рут (в), (54) 
К=1 
Об в = В о», (р)т, (0), (55) 
К=1 


здесь 
1 ; 1 
ор = [Вау ау [8-1] 


о’ (р) = (5-5) а, (р); 
о 
обр) = 7 [4 Круче ()-4,(Р))] 

Характер решений (33)-(50) существенно зависит от типа корней д. Погранслойные первые члены 
этих решений соответствуют краевому эффекту Сен-Венана [4]. В случае мнимых корней д», эти пограничные 
слои имеют слабое затухание. Таким образом, напряженно-деформированное состояние достаточно далеко от 
торцов существенно от них зависит. То есть в этом случае трансверсально-изотропные свойства неоднородного 
материала значительно, по сравнению изотропным материалом сферы, меняют картину напряженно- 
деформированного состояния. В то же время, при действительных или комплексных д, картина напряженно- 
деформированного состояния неоднородной сферы для таких материалов качественно совпадает, различаясь 
скоростью затухания вышеописанных погранслойных решений Сен-Венана неоднородной плиты. 

Из (51) получается, что при удалении от конических сечений 9=0,(/=1,2) решения (33)-(50) 
экспоненциально убывают. 


Поскольку построенные решения удовлетворяют уравнению равновесия и граничным условиям на 
боковой поверхности, вариационный принцип Лагранжа принимает следующий вид [4, 11]: 


Ле - Л (р))би, +(о„» Л. (Ф))би, | м е’?4р =0. (56) 
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Подставляя (52)-(55) в (56) и считая ОЁ, независимыми вариациями, получаем бесконечную систему 


линейных алгебраических уравнений 
Ува =т» (=12,...), (57) 
К=1 

здесь 


1 2 1 
ан = | 0$ (ра (р)е’*ар| Ут, (6, т! (0, ) |+ | о’ (ра (р)е’арх 
1 5=1 =] 


5 (т, ‘био |+ [оь (ру, пруеар (т (@)т, © 


2 


=>. ти (6,)[ 1. (Ра (р)е?ар-+т,(6,)| №, (рда (р)е??ар | 


5=1 


Система (57) всегда разрешима при физически осмысленных условиях, наложенных на правую 
часть (57). Разрешимость и сходимость метода редукции для (57) доказана в [12]. 

Используя малость параметра = , можно построить асимптотические решения системы (57). 

Результаты исследований. Проведен анализ структуры напряженно-деформированного состояния 
радиально-неоднородной трансверсально- изотропной сферы малой толщины при кинематических условиях на 
боковой поверхности. Показано, что, в случае закрепления боковой поверхности, характер решения 
определяется погранслоями. Оказалось, что асимптотическое разложение напряженного состояния начинается с 
решения описывающего краевой эффект Сен-Венана в теории трансверсально-изотропных неоднородных плит. 
В случае трансверсальной изотропии радиально неоднородного материала сферы некоторые погранслойные 
решения затухают весьма слабо, могут проникать глубоко вдали от конических сечений и изменять картину 
напряженно-деформированного состояния. Выведены асимптотические соотношения для перемещений и 
напряжений, позволяющие рассчитать трехмерное напряженно-деформированное состояние радиально- 
неоднородной трансверсально-изотропной сферы малой толщины с любой наперед заданной точностью. 
Показано, что разветвление корней порождает счетное множество новых решений для трансверсально- 
изотропной радиально-неоднородной сферы. 

Обсуждение и заключения. Асимптотический анализ напряженно-деформированного состояния 
неоднородных оболочек, основанный на трехмерных уравнениях теории упругости, позволяет установить 
границы применения приближенных теорий. Выявленный характер поведения решения вдали от торцов для 
разных граничных условий на боковых поверхностях может стать основой для создания уточненных 
прикладных теорий расчета деформирования радиально-неоднородной  трансверсально-изотропной 
сферической оболочки малой толщины. Одним из приложений проведенного асимптотического анализа может 
служить расчет оболочек с тонкими покрытиями, в которых возникает при этом радиальная 
неоднородность [13, 14]. 
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О неоднозначности механической мощности 
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Владимирский электромеханический завод (Владимир, Российская Федерация) 
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Введение. Механические колебания широко распространены в технологических процессах. Приводы машин и 
механизмов преимущественно электромеханические, поэтому механическая реактивная мощность 
трансформируется в электрическую реактивную мощность сети, ухудшая качество электроэнергии. Этим 
обусловлены важность учета механической реактивной мощности и, как следствие, актуальность 
представленной работы. Цель исследования — детализация видов механической мощности при гармонических 
колебаниях. 

Материалы и методы. Изучена литература, в которой освещаются вопросы динамики, кинематики, вибраций, 
преобразования движения в колебательных системах и т. п. Используются теоретические, преимущественно 
математические методы исследования. 

Результаты исследования. Математически осмыслены мощности, развиваемые при упругих деформациях, 
вынужденных гармонических колебаниях инертного тела и колебаниях, связанных с гравитационным 
воздействием, а также реактивная, активная, полная мощности в комплексном представлении и механические 
мощности в векторном представлении. 

Обсуждение и заключения. При механических гармонических колебаниях наряду со знакоположительной 
тепловой мощностью, развиваются знакопеременные реактивные мощности, характеризующие обратимость 
кинетической и потенциальной энергий. Полная механическая мощность удовлетворяет формуле Пифагора. 
Представление о механических реактивных, активной и полной мощностях обобщает соответствующие 
понятия о мощностях из электротехники, и таким образом проявляется электромеханический дуализм. 


Ключевые слова: механическая мощность, кинетическая энергия, потенциальная энергия, комплексное 
представление, векторное представление. 
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Введение. Механическая энергия бывает обратимой (потенциальная и кинетическая), а также 
необратимой (например, тепловая при трении). Временную производную от последней принимают за 
механическую мощность. Отметим, что в силу необратимости тепловой энергии ее производная принимает 
только положительные значения. Вместе с тем производные получают как от потенциальной, так и от 
кинетической энергии. Особый интерес представляют гармонические колебания [1-4], при которых 
производные (мгновенные мощности) будут знакопеременными функциями, что принципиально отличает их от 
тепловой мощности. 

Аналог кинетической энергии в электротехнике — энергия магнитного поля катушки индуктивности, 
аналог потенциальной энергии — энергия электрического поля конденсатора, а аналог механической тепловой 
энергии — тепловая же энергия, рассеиваемая резистором. 

Механические колебания широко распространены в разнообразных технологических процессах [5-8]. 
Приводы машин и механизмов преимущественно электромеханические [9-12], поэтому механическая 
реактивная мощность трансформируется в электрическую реактивную мощность сети, ухудшая качество 
электроэнергии [13]. В этой связи учет механической реактивной мощности имеет немаловажное значение [14], 
и этим обусловлена актуальность представленной работы. 

Материалы и методы. Рассмотрены механические мощности при гармонических колебаниях. В 
качестве литературной базы изучены отечественные и зарубежные источники, в которых освещаются вопросы 
динамики, кинематики, вибраций, преобразования движения в колебательных системах и т. п. Используются 
теоретические (преимущественно математические) методы исследования. 

Результаты исследования 

Мощность, развиваемая при вынужденных гармонических колебаниях инертного тела. Движение 
тела описывается известным выражением: 

х=19шп @Е. 

Соответственно, скорость: 

у=Х = [00$ @{ =И, с0$ @Х. 


Для гармонической величины действующее значение меньше амплитудного в 52: 
| 
ры (1) 


Формула для силы имеет вид: 

` = и = -риау з ах. (2) 
Формула для силы трения: 

и = дх= д16с0$8 @Е. (3) 
Результирующая сила: 


2 . 
Л=Л, + Л, =то’ зтах + шоесозах = 


= 19? +т? о’ 





с05 @# $ш @! 


и та 
и +т2@? Ми’ +т? о? 


Обозначим: 
ф=агои м. (4) 
д 


С учетом этого: 
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= и? +т? в? (созфсоз ах — зп фэт ах) = 1? +т? о? со5(ай +6). 
Очевидно, что 
Е, 2 ре 
= 12| и +779”. 


Действующее значение результирующей силы: 


ыы Гек? + т? 6) 
2 В 


Мгновенная результирующая мощность: 


5= № = и? + та? соз(аи + ф/юсозах = 
= 0, 5? [ит [созф+с0з(2а% + ф)| = 


= РУ [созф+с03(2а#+ф)] = 








= ЕИ (созф- соз2а1 созф-зш 2а/зтф) = 
= ЕТ созф(1+ 05241) — РУ зтфят2ах = р-+а.. (6) 


В электротехнике есть выражение, аналогичное (6), с заменами А — И И - Г. Из него определяют 
активную мощность: 


Р = со$ф. 
Поэтому активную (тепловую) механическую мощность тоже следует определить, как: 
Р = [У со$0. (7) 


Очевидно, что гармонические сила и скорость совершают колебания со сдвигом фаз, равным ф. 


Из вышеназванной формулы электротехники определяют реактивную мощность: 


Р=Шзтф. 
Поэтому реактивную (инерционную) механическую мощность тоже следует определить, как: 
О. = ЕТ зтф. (8) 


Из (6) следует, что под активной мощностью понимается среднее за полпериода значение мгновенной 
мощности, а под реактивной —Щ_ амплитудное значение. В электротехнике аналогично. 
Еще одно обобщение из электротехники — полная механическая мощность: 


$5=ЕТ =. /02+Р2. (9) 


Она примечательна тем, что, с одной стороны, описывается формулой Пифагора, а с другой — равна 
произведению действующих значений гармонических величин. 
Имея в виду (1), (5) и (8), 


пеш? +т? а? 10 та тРах 
О =ЕГзшф= ——= = | (10) 
в Вр: › 








При этом: 
ДХу= та? эт о сов ах = —0,5Р тах зшт 2ах = -Е.У зт2ах = -О, зт24х. (11) 


Это соответствует выражениям (6) и (10). 
Имея в виду (1), (5) и (7), 


и’ +т? о’ 6 и и? 
Р=ЕУ созф = = | (12) 
ий 12 М +т? о? 2 
При этом: 


Ду = Ш со$ в 1@2с0$ ах = 0,5441? 7 (1+с052а4) = ЕТ +с0$2474) = РИ +с0$2а»). (13) 





Это соответствует выражениям (6) и (12). 
Имея в виду (9), (10) и (12), 


ии и’ +т? о’ 1 _ Ро и’ то? 
2 2 | 
Мощность, развиваемая при упругих деформациях. Выражение для силы имеет вид: 
Д. = Ах = М зв аи. (14) 


С учетом (3) результирующая сила равна: 
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Л = Л + Л, = И зах + соз ах = 


"ТТТ К . 
=) К? +7 рек ен 
1+ Г 1? + и’? 
Обозначим: 
ф= ме 
Значит, 


=? + * а (зп фут ах + с08фсоз ах) = [1 + и 08 (и -ф). 


Очевидно, что: 
Е 1? +’. 


Действующее значение результирующей силы равно: 


_ Е, МЕ (15) 


р В 
Мгновенная результирующая мощность: 
5= № = КЕ + и? соз (ай — ф)1осоз ая = 
= 0, 5Ра/? + и? а [соз ф+с05 (2 -— Ф) | = 
= РУ [созф+с0$(2ах—ф)] = 
= ЕУ (созф- соз2а4 созф-+ зш 2ай зтф) = 
= ЕТ созф(1+с0$ 241) + РУ зпфзт2 ах = р-а,. (16) 
Имея в виду (6), (7) и (12), активная механическая мощность равна: 
1? + а? 19 до  _ ига’ 
№2 2 уе + и’ 2 
Принимая во внимание (15), (1), (8) и (16), механическая реактивная (упругая) мощность равна: 
о, гут МН 16 | _ Го. 7) 


№12 №12 ие + и’? 2 


ВР = ЕТ со5ф = 
































При этом: 
Ду= И эт ах оз ах = 0, 5 @2зт 2ах = Е.И зт2ах = О, зт 24. (18) 
Это соответствует выражениям (16) и (17). 
Очевидно, что полная мощность равна: 


$ =АТИ = [02 +Р? Рони уе 
2 
Мощность при колебаниях, связанных с гравитационным воздействием. При отклонении 
подвешенного груза на угол @ возникает момент: 
М = тега. 
Пусть 


@ = а 5 аи. 


@ = 96005 а = бу. й с0$ @1. 


Мгновенная мощность имеет вид: 


Тогда 


4, =Ма = тега, тел соз ай = 0,5та? \1е? эт 29. 
5 0 0 Г 0 


Ее амплитуда и, соответственно, реактивная мощность гравитационного воздействия определяется, как: 


9 = 0, 5тог? 19°. 
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Реактивная, активная и полная мощности в комплексном представлении. В [15] показано, что при 
инертной нагрузке: 


И =Ие"?, 


т 


Мгновенная скорость при этом равна: 


УЕЙ, с05 @Ё = ПАЙ, . 
Формулы ДлЯ действующих значений величин принципиально не отличаются: 
и ре р ро) 
Г=Ие”"", Е=Ее . 
В электротехнике подробно описана особенность комплексного представления: при вычислении 
полной мощности один из перемножаемых векторов должен быть сопряженным. 
э * 
/(л/2 —1 7/2 (л/2+ф-лп/2 Г . ы . 
$ =ЕУ = Ее "Ре"? = рус? = Ус = ЕУ с0зф+ ЕУзшф=Р+/Ю. 
Это выражение для инертной нагрузки. Упругая нагрузка отличается тем, что реактивная мощность 
имеет противоположный знак: 
.* Р ; . . 
5=ЕУ = Ее’ Ву" = рус = Ку № = ЕУ созф- У зтф=Р+ Ю,. 
При этом: 
.* . + 
Р=ВеКУ, О=ВаЕУ. 
Механические мощности в векторном представлении. В основе комплексного представления лежит 


идея вращающихся в комплексной плоскости векторов. Тот же принцип может быть реализован в трехмерном 


Декартовом базисе. 
Из (7)-(9) следует: 


Р=(Е,Г), 9=[Е,У], 5? =(Е,Г) [Е]. 





Математическая абстракция с проекциями вращающихся векторов имеет конкретную материальную 
основу в виде кривошипно-кулисных механизмов. 

Обсуждение и заключения. Математическими методами исследованы мощности: 

— при вынужденных гармонических колебаниях инертного тела, 

— при упругих деформациях, 

— при колебаниях, связанных с гравитационным воздействием, 

— реактивная, активная и полная (в комплексном представлении), 

— механическая (в векторном представлении). 

Показано, что при механических гармонических колебаниях развивается не только 
знакоположительная тепловая мощность, но и знакопеременные реактивные мощности, характеризующие 
обратимость кинетической и потенциальной энергий. 

При этом полная механическая мощность удовлетворяет формуле Пифагора. 

Представление о механических реактивных, активной и полной мощностях является обобщением 
соответствующих понятий о мощностях из электротехники, и таким образом проявляется электромеханический 
дуализм. 
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Введение. Рассматриваются ракета-носитель (РН) в полете и динамические составляющие нагрузок от 
воздействия порыва ветра трапецеидального профиля. Предлагается определить динамические составляющие 

















силовых факторов с помощью аналитических решений для ускорений точек конструкции. Цель работы — 
создание методики выбора продолжительности нормативного порыва, при воздействии которого достигаются 
максимальные нагрузки в сечениях корпуса ракеты-носителя. 

Материалы и методы. Ракета-носитель представляется в виде неравномерной балки. Описание ее колебаний 
приводится к системе независимых обыкновенных дифференциальных уравнений, определяющих движение 
эквивалентной системы осцилляторов. Уравнение колебаний осцилляторов под действием трапецеидальной 
импульсной нагрузки решается методом наложений и сводится к вычислению интеграла Дюамеля. 
Предлагается получать параметры эквивалентной системы осцилляторов по результатам расчета динамических 
характеристик для конечноэлементной модели РН в программе Мазйгап. 

Результаты исследования. Приведены аналитические соотношения для ускорений точек корпуса РН под 
воздействием порыва ветра трапецеидального профиля. Для балочной модели проведены тестовые расчеты 
ускорений по методике, предложенной в настоящей статье. Эти данные сравниваются с результатами 
конечноэлементного моделирования. С помощью аналитических решений построены зависимости, 
определяющие характер изменения величины изгибающего момента для различных сечений ракеты-носителя 
при варьировании продолжительности порыва ветра. 

Обсуждение и заключения. Благодаря представленной методике можно на примере ракеты-носителя строить 
эквивалентную динамическую модель систем с большим числом степеней свободы и получать аналитические 
решения для ускорений точек механической системы при трапецеидальном внешнем воздействии. Данные 
решения применимы для исследования динамических нагрузок. Итоги такого анализа позволяют выбрать 
продолжительность порыва ветра, при которой достигаются максимальные нагрузки в сечениях корпуса РН. 
Вычисления на базе аналитических решений весьма экономичны с точки зрения затраченного времени. Их 
можно использовать в проектных расчетах для предварительной оценки нагружения. 


Ключевые слова: ракета-носитель, балочная модель, осциллятор, корпусные нагрузки, интеграл Дюамеля, 
динамические нагрузки, порыв ветра, аналитические решения, дифференциальное уравнение, колебания, 
изгибающий момент. 
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Введение. При эксплуатации ракеты-носителя (РН) в элементах ее конструкции возникают нагрузки. 
Речь идет о продольных и перерезывающих силах, изгибающих и крутящих моментах. Данные об указанных 
силовых факторах задействуют для прочностного анализа при проектировании новых изделий, 
экспериментальной отработке конструкции [1] и адаптации средств выведения под конкретный пуск [2]. 
Нагрузки подразделяются на квазистатические и динамические. Квазистатические, возникающие в полете, 
вычисляются из условия динамического равновесия РН как твердого тела с учетом допустимых параметров 
траектории выведения. 

Для расчета динамических нагрузок обычно применяются обобщенные балочные модели. Такие 
нагрузки определяются по результатам решения уравнения движения упругой РН, которое в общем случае 
представляет собой уравнение в частных производных. Методы, основанные на использовании разложений 
решений по тонам собственных колебаний конструкции, способны обеспечить высокую скорость вычисления в 
сочетании с достаточной точностью результатов [3]. Как показано в [4], с помощью метода разложения по 
формам собственных колебаний можно перейти к системе независимых обыкновенных дифференциальных 
уравнений. Они описывают: 

— движение центра масс РН, 

— поворот продольной оси РН относительно центра масс, 

— поперечные упругие колебания РН в полете. 

Переход к независимым уравнениям, описывающим упругие колебания РН, означает, что будут 
учитываться распределенные параметры осцилляторов, каждый из которых представляет собой систему с 
одной степенью свободы. Движение такого осциллятора можно рассматривать независимо от других и для 
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каждого из них можно получить решение с помощью известных методов теории колебаний. 

С целью расчета нагрузок в полете рассматривается не вся траектория выведения, а только некоторые 
ее точки, так называемые случаи нагружения, характеризующиеся экстремальным значением отдельных 
параметров, влияющих на нагружение, либо максимальным значением нагрузок на отдельные элементы 
конструкции. Один из наиболее важных случаев нагружения — полет РН в плотных слоях атмосферы [1, 5]. 
Влияние неспокойной атмосферы на нагружение РН можно определить статистическими методами [6, 7] либо в 
рамках консервативного подхода, когда учитываются максимально возможные (с некоторым уровнем 
вероятности) характеристики ветра. В данной статье рассматривается второй подход. В качестве внешнего 
динамического воздействия принимается однократный нормативный порыв ветра. Профиль нормативного 
порыва, характеризующий изменение скорости ветра с течением времени, может быть задан в 
трапецеидальном [8], косинусоидальном или синусоидальном виде [9]. В данной статье рассмотрим движение 
ракеты-носителя под действием трапецеидального порыва. Продолжительность нормативного порыва обычно 
выбирается сопоставимой с периодом низшего поперечного тона колебаний РН'. При этом часто выдвигаются 
требования ее варьирования для достижения максимальных усилий в сечениях РН [8, 10]. Трудоемкость 
расчетов с использованием стандартных программ конечноэлементного (КЭ) анализа обусловлена: 

— необходимостью варьирования параметров внешнего воздействия, 

— большим числом расчетных случаев, 

— множеством вариантов исполнения и конфигураций конструкций на этапе проектных расчетов [9]. 

Цель настоящей статьи — разработка методики выбора продолжительности действия нормативного 
порыва с использованием аналитических решений, полученных для упрощенной динамической модели РН, 
представленной в виде эквивалентной системы осцилляторов. 

Полуаналитический подход с использованием интеграла Дюамеля был успешно применен в [11] для 
проведения гидроупругого анализа судов. В [12] и ряде других работ интеграл Дюамеля используется в рамках 
решения задачи нагружения мостов подвижными нагрузками. В настоящей статье интеграл Дюамеля 
задействовали для аналитических решений реакции ракеты-носителя на кратковременное воздействие порыва 
ветра в полете. 

Материалы и методы. На этапах эскизного проектирования целесообразно использовать плоские 
расчетные схемы для балочных моделей. При простоте и скорости решения они позволяют определить 
параметры движения и внутренние усилия (с приемлемой для этого этапа проектирования точностью) [13]. 
Представим РН в виде упругой балки с переменными по длине массой и жесткостью. Примем обычные для 
сопротивления материалов допущения, в том числе гипотезу о малости упругих деформаций. Для определения 
внутренних силовых факторов в сечении РН используем метод ускорений (перегрузок) [1, 4], который можно 
интерпретировать как метод сечений, адаптированный для динамического расчета. В этом случае внутренние 
усилия находят из условий статического равновесия мысленно отсеченных частей конструкции под действием 
внешних распределенных нагрузок, дополненных силами инерции Д’Аламбера, и искомых внутренних усилий. 
Квазистатические и динамические значения силовых факторов определяются отдельно на основе 
предварительно вычисленных ускорений, а затем суммируются [4]. 

В данной работе рассматривается вопрос определения динамического нагружения РН в поперечном 
направлении под действием порыва ветра, скорость которого направлена перпендикулярно продольной оси РН. 
Принимается, что нагружение в продольном направлении можно рассчитать независимо. В данной работе оно 
не рассматривается. 

С целью определения динамических ускорений РН представляется в виде свободной упругой балки. Ее 
движение изучается в окрестности момента времени, соответствующего рассматриваемому случаю нагружения, 
и описывается в отклонениях от состояния динамического равновесия, в котором находилась РН до ветрового 
порыва, двигаясь по номинальной (невозмущенной) траектории. При этом как бы «замораживаются», 
т. е. принимаются постоянными и равными характеристикам рассматриваемой точки номинальной траектории 
такие параметры, как масса и момент инерции РН, тяга двигателя, угол наклона траектории. Возмущенное 
движение упругой РН исследуется в неподвижной системе координат, связанной с положением, которое РН 
занимала на момент начала расчета. Возмущенное движение будет представлять собой совокупность 
плоскопараллельного движения РН как твердого тела в плоскости, в которой приложена динамическая 
нагрузка, и упругих движений корпуса. Реакцию системы управления не учитываем, т. е. для автомата 
стабилизации предполагаем большое время запаздывания по сравнению со временем приложения 
динамической нагрузки. Вообще, под воздействием порыва ветра вместе с упругими колебаниями корпус РН 
начинает движение как твердое тело. Поперечная составляющая аэродинамической силы, которая при малом 


'Лиходед А. И. Динамика конструкций и определение нагрузок. Королев : Изд-во АО ЦНИИмаш, 2020. 239 с. 
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угле атаки считается ему пропорциональной, меняет значение за счет смещения корпуса в направлении 
действия порыва ветра и поворота относительно вектора набегающего потока. Проекция силы тяжести на 
поперечную ось РН также меняется. Для правильного учета указанных изменений уравнения движения РН 
необходимо интегрировать с уравнениями, описывающими логику работы автомата стабилизации, что 
невозможно на ранних этапах проектирования. Учитывая к тому же значительную массу и момент инерции РН, 
будем считать малыми: 

— угол поворота РН как твердого тела за время расчета, 

— скорость смещения центра масс РН в направлении действия порыва ветра. 

Это позволяет не учитывать влияние указанных выше изменений. Приращение угла атаки (а 
следовательно, и поперечной составляющей аэродинамической силы) считается зависящим только от величины 
скорости порыва ветра, заданной в виде функции времени. Таким образом, с учетом принятых допущений 
аэродинамическая нагрузка в поперечном направлении представляет собой распределенную по длине балки 
нагрузку с коэффициентом пропорциональности, зависящим от времени. Закон распределения 
аэродинамической нагрузки по длине РН определяется экспериментально и считается заранее известным. Закон 
изменения коэффициента пропорциональности (угла атаки) от времени определяется выбором профиля 
нормативного порыва ветра. 

Движение РН, смоделированной в виде упругой балки, можно описать, используя известное уравнение 
вынужденных поперечных колебаний балки, записанное с учетом гипотезы Фойгта: 


92 д д2 
то + (1+2) 2. [Вод] = ч.5, (1) 
где т(х) — погонная масса; В(х) — изгибная жесткость; 4 (х, ь — распределенная внешняя нагрузка; й — 
коэффициент трения. 


Данное уравнение должно быть дополнено граничным условием: внутренние усилия в начальном и 
конечном сечении равны ьы. Значит: 


> [В(х о =0, В(х уе — 0, прих = 0,х = [, (2) 


где [ — длинаРН. 
В качестве внешней распределенной нагрузки в данной работе принимается аэродинамическая сила, 
которую можно представить в виде произведения функций: 

4,0 = ВУ), (3) 
где Е (Е) — функция, определяющая временную изменчивость аэродинамической силы и изменяющаяся по 
трапецеидальному закону в соответствии с принятой в настоящей работе моделью порыва ветра; У„(х) — 
функция распределения аэродинамической силы по длине РН. 

Рассмотрим свободные колебания РН без учета сил трения (при 49(х, Е) = 0, В = 0). Подставим метод 
разделения переменных у(х, Е) = /(х)`4(®). В этом случае от уравнения (1) с граничными условиями (2) 
можно прийти к обыкновенному дифференциальному уравнению 











42 ры 
ча [ВР]  р2т(х) ©) = 0, (4) 
с краевыми условиями: 
< [ВС Е = 0, В(х О =, прих = 0х = [. (5) 


Решение (4) с а (5) представляет собой классическую задачу Штурма — Лиувилля. Решая ее, 
можно найти набор собственных форм Д(х) и собственных частот р; рассматриваемой балки (] = 1,2,...). 
Известно”, что некоторые решения системы (4) соответствуют нулевым собственным частотам. Формы, 
соответствующие нулевым собственным частотам, определяют поступательное движение РН как твердого тела 
вместе с центром масс и вращение вокруг центра масс: [_. = 1, К =Х — хс, где хс — координата центра масс 
РН. 

Следует отметить, что массовые и жесткостные характеристики РН чаще всего имеют кусочно- 
постоянный характер распределения. В этом случае уравнения вида (1) и (4) должны быть записаны отдельно 
для каждого однородного балочного участка с граничными условиями на стыках участков, как при выводе 
соотношений метода начальных параметров [1, 14]. В данной работе указанная запись опущена, т. к. расчет 
динамических характеристик (модальный анализ) конструкции проводится численно, с использованием метода 
конечных элементов. 

Представим вынужденные колебания упругой балки, моделирующей корпус РН, в виде разложения по 
формам собственных колебаний. Примем, что ось жесткости балки проходит через ее центр масс. Для 
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перемещения точек оси РН запишем: 

Убе, = ус(®) + 90 (х — хе) +.) а; ©, (6) 
где ус(Е) — перемещение центра масс балки; 9(Е) — угол поворота оси недеформированной балки; [, — форма 
собственных колебаний балки, соответствующая тону с номером }; 4;(Е) — обобщенная координата, 


соответствующая тону с номером ]; М — число учитываемых упругих тонов. 
После подстановки (6) в (1) и применения процедуры Бубнова — Галеркина можно прийти к системе 
обыкновенных дифференциальных уравнений с постоянными коэффициентами: 


тус(®) =0,, 
196) = 05, (7) 
т; (4; (©) + 2п;4;(® +274;(0) = 90 Ч =Ъ2,...,№ 
Здесь т — масса РН; [ — момент инерции относительно оси, проходящей через центр масс РН 
перпендикулярно плоскости вращения; 7; — приведенная (обобщенная) масса для -го тона колебаний и 


определяемая формулой: т; = т т(х) }? (х) ах. 
Обобщенные силы в выражении (7) с учетом (3) определяются следующим образом: 


7Й 
90, == В(0 у У, 0) ах = ВОО, (8) 
7Й 
05 = В(® Л Ух) (х — ходах = В(®Ма, (9) 
7Й 
90; =В@) т У„ (5) фах = В(@) 0. (10) 
Здесь О а —Щ_ максимальное значение главного вектора поперечной аэродинамической нагрузки; Ма —— 


максимальное значение главного момента поперечной аэродинамической нагрузки, приведенного к центру масс 
РН; 90; — максимальное значение обобщенной силы, соответствующей обобщенной координате 4). 

Первые два уравнения в (7) определяют закон изменения ускорений точек РН в процессе 
поступательного и вращательного движения РН как твердого тела. Последнее уравнение в (7) определяет закон 
движения эквивалентной системы осцилляторов. 

Рассмотрим движение одного осциллятора под действием трапецеидальной внешней нагрузки, которая: 

— возрастает от нуля до Оу за время 6, 

— сохраняет постоянное значение за время 6, 

— падает до нуля за время 6. 

Для удобства опустим индексы, характеризующие номер тона. Тогда дифференциальное уравнение 
движения осциллятора с учетом (10) перепишем в виде: 


4+2п4 + р?4 = =%5(0, (11) 
где т — масса осциллятора; р — угловая частота собственных колебаний, выраженная в радианах в секунду; и 
— коэффициент затухания (определяет демпфирующие осциллятора); 0 (Е) = О5Е(®) — закон изменения 
внешней нагрузки. 


Представим функцию В (Е) как совокупность четырех линейных функций: 


В. (© =б *%, В2( =-6 (Е - 6), 
В-(® =-6 1(Е- 6-0), В.(© = 6 "(Е- 26-0). 
Соответственно, внешняя нагрузка О(Р) — это совокупность четырех линейных нагрузок (0;:(Е) = 


ОВ; (©, (Е = 1,2,3,4). Нагрузка 0.(®) прикладывается с момента & = 0; 0>(Е) — с момента & = 5; 03(® — с 
момента & = 0 +5; О.(Е) — с момента Е = 0 + 25. 

Для определения реакции системы на внешнее воздействие разобьем все время действия нагрузки на 
четыре интервала (рис. 1). 
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Рис. 1. Схема внешней трапецеидальной нагрузки 
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В соответствии с методом наложений? представим в виде интеграла Дюамеля реакцию 
рассматриваемой линейной системы на внешнее воздействие как сумму реакций на совокупность независимо 


приложенных элементарных импульсов: 
Е 


а(® = | Р(Е— 8) (8) а, 


0 
здесь Р(Е — 9) — смещенный временной закон изменения внешнего воздействия, а У (9) характеризует реакцию 


системы на единичное импульсное воздействие. Реакцию механической системы 41(Р) на линейно 
возрастающую нагрузку Р(Е) = КЕ можно выразить через реакцию системы на единичный импульс У(9) и на 


внезапно приложенную единичную нагрузку У' (9): 
р с 


О | у. (8), У®= | у(з)аз. 


0 0 
Реакция на единичный импульс для механической системы с одной степенью свободы, имеющей 


коэффициент затухания п, будет иметь вид": 


1 
У (9) = 
инет 


Введем обозначение для частоты затухающих колебаний р, = \/р? — п?. Вычислим реакцию системы 


= та ( р? — п25). 


на внезапно приложенную единичную нагрузку: 
Е 


1 п 
У, (9) = [уфа =—- Ё ет (со5Ф + — 5шФф 5) 
1 ) тра 1 р, г 
Реакция на линейно возрастающую нагрузку: 
а(® = |, (9)а9 = — [-1+е-"* соз(р1#)] +2 (= _ 1) е-"* зт(р: В) + :} (12) 
1 071 тр? (р? Ра и т 1 . 
Примем обозначения: Л. = Ё -— 6, Д›= Д, — 0, А. =Д, — 6. 


Введем функцию, содержащую гармонические члены решения (12): 
п п 1 
4(9) = 2 с (соз(ра 9) + — зт(ру :9)) е-"8 — — зт(р! 9) е-"8. 
р Рл Рл 


Обозначим 4. перемещение под действием статически приложенной к системе силы 0% = 4.тр? и 
учтем, что К = 05/6. Суммарное перемещение осциллятора под действием комбинации нагрузок 0;(Е) на 
каждом временном участке будет суммой соответствующих решений (12). Реакция системы, описываемой 
уравнением (11), на внешнее воздействие трапецеидального профиля будет иметь вид: 


Чстб ЧА(® - 2п/р? +8 при 0 < Е< 65, 
а(® = 46 Ча(®) - а(4,) + 6] приб < Е < 0+5, (3) 
аб Ча(е) - а(4,) - а(4,) + 2.п/р? -4.] приб +6 < Е< 0+ 256, 
Чеб-(а(®) — а(№) — а(Аь) + а (дз) В 
Закон изменения ускорений осциллятора может быть получен двойным дифференцированием по 
времени выражения (13). Введем функцию 9(9) = к (а (9). Дифференцируем и запишем результат с учетом 


номера тона колебаний: 
п; п; 1 
9:09) = [сз — 22) (2 т (обр -9) + —^ зт(рь ; :9)) — —эт(ра; :9)) _ 
р; Р1 Р1 


—2п; (2 ы (-ра; зт(р; .9) + п; соз(р ; .9)) — соз(р: ; :8)) е-—й8. (14) 


Перепишем выражение для частоты затухающих собственных колебаний: 


раз = |2? п. и 


Для определения ускорений точек РН дважды продифференцируем (6): 
У(Е,х) = Ус (®) +9(0(х — х) +. 60а. (16) 
Выразим ускорения обобщенных координат из первых двух уравнений системы (7) с учетом (8) и (9). 
Представим ускорения точек оси РН, которые она приобретает, двигаясь как твердое тело: 


а(®) = ус (© +90 (х — хе) = В(® (@а/т + Мо/П. (17) 


- Бидерман В. Л. Теория механических колебаний. М. : 0В$$, 2017. 416 с. 
Яблонский А. А. Курс теоретической механики. М. : Интеграл-Пресс. 2007. 603 с. 
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Кроме того, учтем закон изменения функции В (Е) от времени: 


5-м при 0 < Е < 5, 
в(0 = 1 при б < Е < 9-65, 
—6 (Е- 26-0) при 0 +6 << 0+ 25, 
0 при > 0 + 26. 
Получим: 
9 1(Оа/т + Ма/Г) при 0 < Е< 65, 
а(г) = Оа/т + Ма/1 приб < Е < 9 +65, (18) 
—5 Е - 26 -0)(0./т+М./Г) при0+6 < Е<0+26, 
0 приё > 0 + 26. 


Для ускорений обобщенных координат, соответствующих упругим тонам колебаний, в результате 
двойного дифференцирования выражения (13) с учетом (14) получим: 
5 ‘Чу 9:0) при 0 < &< 65, 


О (ях — 9/(А,)) приб << 0+5, 
5-е (о — 9 (А) — 9, (45) при +5 << 0+ 25, 


5-4) (9/0 — 9) — 9,45) + 9,(43)) приё>0+25. 

С учетом (16) и (17) закон изменения ускорений точек упругой оси РН под действием порыва ветра 
трапецеидального профиля будет иметь вид: 

Ув, х) =а(® +1. 2 0)4,(0, (20) 
где а(2) определяется выражением (18), 4;(Е) — выражением (19). 

При исследовании упругих колебаний РН демпфирование традиционно принимается на основании 
данных, полученных по результатам натурных динамических испытаний и представленных в виде значений 
логарифмических декрементов 0;. Тогда коэффициент, определяющий параметр демпфирования и входящий в 
выражения (14) и (15), можно вычислить по формуле: 

р, 

Входящее в формулу (19) изменение обобщенной координаты 4; под действием статически 
приложенной обобщенной силы @о; определяется выражением: 

1 00) __ 1 р абледах 
рт) р) утодах 

Принимая во внимание (14), (15), (18), (19), (21), (22), формула (20) — это аналитическое выражение, 
определяющее функции изменения ускорений по времени для точек оси РН под действием внешней 
поперечной аэродинамической силы, изменяющейся по трапецеидальному закону с учетом влияния 
диссипативных сил. 

Зная закон изменения ускорений, можно по известным методикам [4] определить динамические, а 
затем и суммарные корпусные нагрузки, действующие в сечениях РН. Вычислим изгибающий момент, 
обусловленный инерционными силами, возникающими за счет упругих колебаний корпуса РН: 


4, (#) = (19) 


Чст Я. — (22) 


М =-У-. Мда. (23) 

Здесь М;(х) — функция распределения по длине РН единичного (при ускорении 4;(Е), равном единице) 
изгибающего момента для -го тона колебаний. Ее можно найти по формуле: 

М; (х) = Г тор (х)ах ах. (24) 


Как отмечалось выше, для проведения модального анализа в данной работе использовался метод КЭ. 
Данный подход обусловлен тем, что на практике у динамической модели РН достаточно сложная структура. 
Она включает подконструкции, и нельзя пренебрегать их собственной динамикой. Подконструкции могут 
крепиться к корпусу РН в одном сечении либо располагаться параллельно продольной оси РН и иметь 
несколько точек крепления. В этом случае расчет динамических характеристик в континуальной постановке 
является сложной математической задачей. Кроме того, динамические модели отдельных подконструкций 
представляются предприятиями-разработчиками в конденсированном (матричном) виде в формате Мазгап. По 
этой причине для расчета динамических характеристик конструкции будет оптимальным задействовать пакет 
программ инженерного анализа Мазгап. 

Однако применение стандартных программ КЭ-анализа для расчета динамического нагружения 
сопряжено с определенными трудностями. К ним относится необходимость предварительно строить 
эквивалентную модель внешней аэродинамической нагрузки, пригодную для использования в программе КЭ- 
анализа [15], что само по себе достаточно сложно. Кроме того, есть трудности, связанные с обработкой 
результатов расчета. Применение для анализа результатов функционала постпроцессора чрезвычайно 
трудоемко, требует большого числа ручных операций. Другой путь предполагает применение дополнительного 
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программного обеспечения для обработки массивов числовых данных большой размерности. В данной работе 
предлагается подход, в рамках которого программа КЭ-анализа используется только для модального анализа. 
Динамическое нагружение в этом случае рассчитывается с помощью специально разработанного программного 
обеспечения, позволяющего варьировать внешние нагрузки и автоматически обрабатывать результаты расчета. 

Для сложной конструкции РН уравнения (7) сохраняют свой вид [3]. Стандартная информация вывода 
программы Мазфап может быть базой для получения параметров эквивалентной системы осцилляторов, массы 
и момента инерции РН, а также для вычисления обобщенных сил, входящих в третье уравнение системы (7). 
Для формирования левой части третьего уравнения в (7) необходимы значения собственных частот р; (КаЧ1апз) 
и обобщенных масс т; (Сепегайте4 таз$). Для определения обобщенных сил в правой части третьего 
уравнения в (7) необходимы функции собственных форм колебаний Д.(х) (Е1вепуесюг). При вычислении 
динамических инерционных нагрузок по (23) вместо (24) удобнее воспользоваться единичными инерционными 
нагрузками (силами и моментами), которые выводятся программой Мазгап после стандартного запроса усилий 
при расчете собственных форм и частот. Единичные инерционные нагрузки выводятся отдельно для каждого 
тона колебаний, причем помноженные на квадрат собственной частоты, что следует учитывать для их 
корректного использования. 

Результаты исследования. Для проведения тестовых расчетов рассмотрена ракета-носитель среднего 
класса тандемной схемы. Вид функций погонной изгибной жесткости В(х) и массы т(х), а также 
распределения сосредоточенных масс тр(х) = Теоср,А(Х —х,) по длине рассматриваемой РН представлен 


на рис. 2 (через А обозначена дельта-функция Дирака). 


т(х), кг/м 


Тсоср (%), кг 


кг.Н 


во, Ета 


в) 


Рис. 2. Распределение массовых и жесткостных характеристик по длине РН: 
погонная масса (а), сосредоточенная масса (6), изгибная жесткость (в) 
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Для проверки полученных аналитических решений при воздействии трапецеидальной внешней 
нагрузки рассчитаны динамические ускорения РН по формуле (20) с учетом (14), (15), (18), (19), (21), (22). При 
этом учтены 5 упругих тонов собственных поперечных колебаний РН. Характеристики р;, т; и //(х) были 
получены по результатам расчета в программном комплексе М$ЗС Мазфап с использованием 
последовательности решения для модального анализа собственных колебаний (ЗОГ 103). Динамическая 
конечно-элементная модель РН представлена в виде набора балочных элементов с различными инерционными 
и жесткостными характеристиками. К балочным упруго или жестко присоединены: 

— элементы, описывающие инерционные свойства приборов, агрегатов, частей конструкции блоков; 

— конденсированные модели отдельных блоков, представленных в цифровом матричном виде. 

Конечноэлементная модель РН включает около 1000 составляющих. Элементы, моделирующие 
закрепления, не использовались для сохранения способности движения РН как твердого тела. 

По той же конечноэлементной модели РН в программном комплексе М$С Мазап рассчитаны 
ускорения с использованием последовательности решения БОГ 119, применяемой для модального анализа 
переходных процессов. При этом в модальном разложении были учтены все тона собственных колебаний в 
диапазоне до 100 Гц. Аэродинамическая нагрузка представлена поперечными погонными нагрузками, 
распределенными по всем балочным элементам, моделирующим корпус РН. Кроме того, задан 
трапецеидальный закон изменения аэродинамической нагрузки по времени. 

На рис. 3 представлены результаты сравнительного анализа полученных с помощью двух различных 
подходов ускорений ао(#) = У(Е, хо) некоторой точки оси РН с координатой х = хо. Видно, что хорошо 
согласуются два решения: 

—щ полученное для ускорений по упрощенной модели РН на основе аналитических соотношений, 
приведенных в настоящей работе; 

— численное, полученное по полной конечноэлементной модели РН. 


ад(®, 30,0 
м/<° 20,0 
10,0 
0,0 - 
-10,0 
-20,0 
-30,0 





0,0 0,5 1,0 1,5 2,0 


Рис. 3. Сравнение поперечных ускорений, полученных с помощью двух различных подходов: ——— — конечноэлементное 
моделирование, ® ® ® ®® — аналитическое решение для упрощенной модели 


Выбранное из временного процесса пиковое значение изгибающего момента определяет уровень 
эквивалентных усилий, принимаемых для проведения прочностного расчета, и выступает изменяющимся 
параметром при варьировании параметров внешнего воздействия в поперечном направлении [16]. 

На основе значений динамических ускорений получены зависимости (рис. 4). Они показывают, как 
изменение пиковых значений изгибающего момента М для различных сечений РН связано с параметром 6, 
характеризующим продолжительность действия порыва. На рис.4 значения изгибающего момента 
представлены в виде безразмерных величин М*. Они рассчитаны делением размерного изгибающего момента 
на максимальное значение для данного сечения (например для сечения х = 0,3Ё максимальное значение 
Мах = 4,6 : 105 Н-м), найденное при варьировании величины 0 во всем рассматриваемом диапазоне. 
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Рис. 4. Зависимость значения безразмерного изгибающего момента М* от продолжительности действия порыва ветра 0*, 
выраженной в долях периода Т; для различных сечений РН: зона 1 (х = 0,0 — 0,2/Г,) (а); зона 2 (х = 0,2[ — 0,41) (6); зона 3 
(х = 0,4[ — 0,551) (в); зона 4 (х = 0,55 — 0,75Г.) (г); зона 5 (х = 0,75[ — 0,91) (9); зона 6 (х = 0,9[ — Г) (е) 


Итак, продолжительность действия порыва на графиках представлена безразмерной величиной 0*, 
полученной путем деления параметра 9 на период Т, первого тона колебаний РН. Каждая линия графика 
соответствует одному сечению РН. Все сечения РН сгруппированы по характеру функции М*(0*) и приведены 
на различных графиках, а длина РН соответствующим образом разделена на зоны (рис. 5). 
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Рис. 5. Зоны сечений РН 
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Из рис. 4 видно, что для первой зоны максимум изгибающего момента достигается уже при значениях 
0 не менее 15 % от периода Т, первого тона колебаний. Для второй и третьей зон — несколько больше 
половины периода Т,. Дальнейшее увеличение значений 0 не влияет на величины максимальных значений 
изгибающего момента. Для зон 4—6 максимум значения изгибающего момента оказывается локальным и 
располагается в зоне значений 09, близких к величине половины периода Т, первого тона колебаний РН. 
Полученные результаты полностью соответствуют итогам конечноэлементного моделирования, проведенного 
ранее°. 

Таким образом, для получения максимальных значений изгибающего момента в сечениях 
рассматриваемой РН необходим динамический анализ поведения конструкции с внешним воздействием в виде 
порыва ветра, продолжительность которого определяется параметром 0, близким по значению к половине 
периода первого тона колебаний РН. 

Обсуждение и заключения. С помощью метода наложений получены аналитические решения, 
описывающие движение системы с одной степенью свободы, на которую влияют сила трения и внешняя сила, 
меняющаяся по трапецеидальному закону. Приведена методика применения полученных аналитических 
решений для систем с большим числом степеней свободы. Показано хорошее совпадение двух видов решений: 

— аналитических, для ускорений точек РН, найденных по упрощенной модели; 

— численных, полученных по полной конечноэлементной модели РН. 

Показано, что аналитические решения могут применяться для анализа динамических силовых факторов 
с целью выбора продолжительности порыва ветра, при воздействии которого достигаются максимальные 
нагрузки в сечениях корпуса РН. Аналогично можно анализировать перегрузки, которые достигаются в 
сечениях РН (например, в точках установки систем измерения). 

Кроме того, на основании предложенной методики можно построить полный цикл предварительного 
расчета нагрузок в случае, когда возможно аналитическое представление внешней динамической нагрузки. 
Анализ нагружения на основе аналитических решений весьма экономичен с точки зрения расчетного времени и 
может быть прекрасной альтернативой конечноэлементному моделированию на этапе проектных расчетов, 
когда исследуется большое число вариантов комбинаций внешних нагрузок и конфигураций разрабатываемой 
конструкции. Конечноэлементный анализ подробной модели в этом случае можно использовать как 
уточняющий итоговый расчет. 
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Введение. Для поиска эффективных методов отделочно-зачистной обработки длинномерных деталей 
целесообразно создать роторно-винтовую технологическую систему проходного типа. Основной ее рабочий 
элемент — винтовой ротор, представляющий собой сочетание плоских элементов различной формы и 
типоразмеров, разнонаправленных по отношению к винтовым линиям по периметру. 
Цель настоящего исследования — обоснование интенсивности процесса обработки в устройствах, оснащенных 
винтовым ротором. 
Материалы и методы. В качестве основного параметра, определяющего интенсивность процесса обработки 
деталей в винтовых роторах, принят съем металла. В рамках представленного исследования обработка 
выполнялась на экспериментальной роторно-винтовой установке. Обрабатывающая среда состояла из 
формованных абразивных гранул марки ПТ 10х10. Условия исследования: 

— объем загрузки (без детали) — 60 %, 

— скорость вращения ротора — 50 об/мин, 

— время обработки — 30, 60, 90 мин, 

— углы наклона оси винтового ротора — 0° и 5°. 
Влияние режимов и условий обработки на интенсивность процесса рассматривалось на пластинчатых образцах 
с размерами 80х10х1 мм из алюминиевого сплава Д16Т. Для определения съема металла образцы до и после 
обработки взвешивали на аналитических весах ОВайз АХ223. 
Результаты исследования. Представлены закономерности съема металла с образцов при различных 
продолжительности обработки, расположении образцов в рабочей зоне винтового ротора с разными углами 
наклона и скорости его вращения. 
Обсуждение и заключения. Установленные в ходе исследований закономерности свидетельствуют об 
эффективности роторно-винтовых технологически систем для решения задач, связанных с отделочной 
обработкой длинномерных деталей. Основные факторы, обеспечивающие управление такой обработкой в 
устройствах с винтовым ротором и влияющие на его интенсивность (и, как следствие, — на 


производительность процесса): скорость вращения винтового ротора, форма периметра и угол наклона оси. 


Ключевые слова: отделочно-зачистная обработка, винтовой ротор, длинномерные детали, обрабатывающая 


среда, интенсивность процесса, съем металла. 
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Введение. Поиск эффективных путей отделочно-зачистной обработки длинномерных деталей [1-6] 
предопределил создание для решения этой задачи роторно-винтовой технологической системы проходного 
1,2 
типа (рис. 1) [7-15] >^. 






























































Я 
57Р | 
52 (. й 
7 Е, 
Выход ТХ 
| 
3 Т 7 т 2 


Рис. 1. Технологическая схема отделочно-зачистной обработки длинномерных деталей в винтовом роторе: 
1 — установочная рама, 2 — проушины, 3 — регулировочные винты, 4 — двигатель, 5 — бак-отстойник, 6 — приводные 
валы, 7 — роликовые опоры, 8 — обечайки, 9 — ролики, 10 — ободы, 11 — винтовой ротор, 12 — прокладка, 
13 — прокладка с прорезью, 14 — система трубопроводов, 15 — деталь, 16 — обрабатывающая среда, 17 — контейнер 


? Устройство для вибрационной обработки длинномерных деталей : патент 2228252 Рос. Федерация : В24АВ 31/06/ А. П. Бабичев, 
И. А. Бабичев, Г. В. Серга. № 2002135225/02 ; заявл. 25.12.2002 ; опубл. 10.05.2004, Бюл. № 13. 7 с. 

ь Агрегат для отделочно-зачистной обработки длинномерных деталей : патент 2750922 Рос. Федерация : В 24В 31/06 / В. А. Лебедев, 
Л. М.Р. Аль-Обайди, Н. С. Коваль, Г. В. Серга. № 2020119703 ; заявл. 15.06.2020 ; опубл. 06.07.2021, Бюл. № 19. 27 с. 
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Основным рабочим элементом этого типа технологической системы является винтовой ротор. Его 
периметр выполнен в виде сочетания плоских элементов различной формы и типоразмеров, разнонаправленных 
по отношению к винтовым линиям по их наружной поверхности. Обрабатывающая гранулированная среда 
находится внутри вращающегося винтового ротора. Увлекаемая плоскими элементами, в определенный момент 
под действием силы тяжести она лавинообразно скатывается вниз. При этом верхние слои масс загрузки 
вращаются вокруг собственных осей, втягивая близлежащие слои частиц масс загрузки. В результате частицы 
масс загрузки перекатываются относительно друг друга. К тому же при скатывании лавиной по неровной 
поверхности нижележащих слоев происходят удары небольшой силы и скольжение частиц масс загрузки. 
Таким образом, деталь обрабатывается истиранием, царапаньем и в небольшой мере ударами. Технологический 
эффект обработки зависит от габаритов винтового ротора, масс частиц обрабатывающей среды, коэффициента 
заполнения рабочей камеры, скорости вращения, геометрии винтового ротора и ряда других факторов. 

Цель представленного исследования — обоснование интенсивности процесса обработки в устройствах, 


оснащенных винтовым ротором. 
Материалы и методы. В качестве основного параметра, определяющего интенсивность процесса 
обработки деталей в винтовых роторах, принят съем металла. 
Для обработки использовали экспериментальную роторно-винтовую установку, изготовленную на базе 
токарно-винторезного станка модели 1К625 (рис. 2). 
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Рис. 2. Общий вид экспериментальной роторно-винтовой установки 


Обрабатывающая среда состояла из формованных абразивных гранул марки ПТ 10х10. Условия 
исследования: 

— объем загрузки (без детали) — 60 %, 

— скорость вращения ротора — 50 об/мин, 

— время обработки — 30, 60, 90 мин, 

— углы наклона оси винтового ротора — 0° и 5°. 

Влияние режимов и условий обработки на интенсивность процесса рассматривалось на пластинчатых 
образцах с размерами 80х10х1 мм из алюминиевого сплава Д16Т. Из этого материала выполнены заготовки в 
виде длинномерного проката, требующего отделки и зачистки перед запуском в производство. 

В рабочей зоне экспериментального винтового ротора по длине рабочей зоны определялись три 
сечения А, В, С (рис. 3), удаленные от левого торца винтового ротора на расстояния 1, [4 [>. 
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Рис. 3. Схема установки образцов на оправке в роторе: 
1 — образцы; 2 — оправка; 3 — резцедержатели суппорта токарного станка 
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Образцы закреплялись на длинномерном четырехгранном полом стержне и устанавливались по схеме, 
представленной на рис. 3. Расстояния к [ [5 соответствуют положению сечений рабочей зоны ротора. Для 
оценки характера воздействия обрабатывающей среды на различные поверхности оправки в каждом положении 
(сечении) обрабатывалось по 4 образца. Их закрепляли на гранях оправки таким образом, что образцы [и 2 
находились в зоне активного воздействия обрабатывающей среды (со скользящим слоем), а образцы 3 и4 — в 
пассивной зоне (рис. 4). 
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Рис. 4. Схема расположения образцов на оправке в сечении: 
1,2, 3, 4 — образцы; 5 — винтовой ротор; 6 — обрабатывающая среда 


Оправка с образцами устанавливалась в резцедержателе суппорта токарного станка и вводилась в 
рабочую зону винтового ротора для обработки. 

Для определения массового съема металла образцы взвешивали до и после обработки на аналитических 
весах Овал АХ223. Данная модель предназначена для статических измерений массы веществ и материалов с 
точностью до 0,001 г. Съем металла контролировали после тщательной мойки и сушки образцов. 

Удельный съем металла с единичной площади поверхности, равной площади квадрата упаковки гранул 


обрабатывающей среды, определялся из соотношения 
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где ДО — съем металла с поверхности, г; & — время обработки, с; 5 — площадь обрабатываемой поверхности 
образцов, мм?; 5, — площадь квадрата упаковки гранул обрабатывающей среды, мм”. 

Продолжительность обработки фиксировалась с помощью реле времени (таймер) модели ВЛ-45УХЛ41. 

Результаты экспериментов определялись как среднеарифметическая разность соответствующих 
значений измерений образцов или партии образцов до и после обработки. 

Результаты исследования. На рис. 5—7 приведены результаты экспериментальных исследований 
интенсивности процесса обработки в винтовом роторе. 
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Рис. 5. Изменение съема материала с образцов Д16Т в зависимости от продолжительности обработки и расположения 
образцов на оправке: с углом наклона 0° (а); с углом наклона 5° (6) 
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Рис. 6. Изменение съема металла с образцов Д16Т в зависимости от их расположения 
в рабочей зоне винтового ротора (за 90 мин): с углом наклона 0° (а); с углом наклона 5°(6) 
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Рис. 7. Изменение удельного съема металла с образцов Д16Т в зависимости от их расположения на оправке (за 90 мин): с 
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Рис. 8. Изменение съема металла с образцов Д16Т в зависимости от загрузки рабочей зоны винтового ротора 


обрабатывающей средой и скорости вращения (за 30 мин) 
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Обсуждение и заключения. Результаты исследований позволяют сделать ряд утверждений. 

1. Вид кривой, описывающей изменение величины съема металла с поверхности от продолжительности 
обработки, близок к линейному (рис. 5). Это не зависит от расположения образцов в рабочей зоне винтового 
ротора, его конструктивного исполнения и угла наклона по отношению к оси вращения. 

2. Кинетика движения частиц и, как следствие, интенсивность съема металла в различных сечениях 
рабочей зоны винтового ротора определяются углом наклона пластин, создающих форму рабочих секций 
ротора. При их одинаковой форме (рис. 3) обеспечивается практически одинаковая интенсивность обработки. 
Разница величин съема в трех сечениях составила 3 %. 

3. Интенсивность обработки (рис. 5—7) зависит от положения поверхности образцов по отношению к 
направлению движения гранул скользящего слоя (рис. 4). Сравнивались два случая съема металла. Первый: с 
образцов 1 и 2, установленных в активной зоне воздействия обрабатывающей среды. Второй: с образцов 3 и 4, 
установленных в пассивной зоне, характеризующейся слабым энергосиловым взаимодействием с 
обрабатываемой поверхностью. В первом случае съем металла вдвое выше. 

4. Изменение угла наклона оси винтового ротора по отношению к оси его вращения увеличивает 
интенсивность процесса обработки образцов в среднем сечении (В-В) в 1,2 раза и снижает ее в рабочих зонах 
крайних секций (рис. 5—7). Это обусловлено изменением кинематики движения обрабатывающей среды в 
рабочей зоне винтового ротора. При каждом полном обороте обрабатывающая среда стремится переместиться в 
среднюю часть рабочей зоны ротора (сечение В-В). Так обеспечиваются: 

— постоянство давления объема среды на образцы; 

—щ активация ударно-импульсного воздействия гранул обрабатывающей среды на поверхность 
образцов. 

Кроме того, при изменении угла наклона оси винтового ротора снижается осевое давление 
обрабатывающей среды на торцевые стенки. Это очень важно для установок проходного типа, т. к. исключает 
высыпание гранул из отверстий, предназначенных для прохождения длинномерных деталей через рабочую 
зону. 

5. Скорость вращения винтового ротора — важный технологический фактор, определяющий 
интенсивность процесса обработки в роторно-винтовой установке. Из рис. 8 видно, что при различной степени 
загрузки величина съема материала со всех образцов растет с увеличением скорости вращения винтового 
ротора до 50 об/мин. При повышении скорости до 70 об/мин показатель снижается. Это вызвано усилением 
центробежных сил. Их действие приводит к тому, что некоторый объем частиц обрабатывающей среды 
постоянно находится у стенок ротора. В результате сокращается количество скользящих слоев 
обрабатывающей среды, и, следовательно, снижается интенсивность воздействия на поверхность образцов. 

Установленные в ходе исследований закономерности и технологические возможности процесса 
обработки в винтовом роторе свидетельствуют о приемлемости и эффективности применения роторно- 
винтовых технологических систем для решения задач, связанных с отделочной обработкой длинномерных 
деталей. 

В проведенных изысканиях основными рассматриваемыми факторами были: скорость вращения 
винтового ротора, угол наклона его оси и форма периметра. Они обеспечивают управление отделочной 
обработкой в устройствах с винтовым ротором, влияют на его интенсивность, а значит, и на 
производительность процесса. 
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Введение. Работа посвящена исследованию триботехнических характеристик экспериментальных твёрдых 
сплавов с модифицированной связкой при трении без смазки по труднообрабатываемым в процессе резания 
материалам — нержавеющей стали и титановому сплаву. Целью исследований является оценка процесса 
фрикционного взаимодействия для каждой пары трения по ряду параметров и определение на основании 
установленных трибопоказателей оптимальных сочетаний «экспериментальный твёрдый сплав — 
конструкционный материал». 

Материалы и методы. Трибологические испытания твёрдых сплавов проводились по схеме трения «цилиндр- 
диск» для различных скоростей скольжения и температур при постоянной нагрузке без применения смазочного 
материала. Сравнение процесса фрикционного взаимодействия производилось по силе трения, объёмному 
износу и шероховатости дорожек трения на контртеле. В качестве материалов контртел применялись 
нержавеющая сталь 12Х18Н9Т и титановый сплав ВТЗ-1. Определение устойчивости экспериментальных 
составов к абразивному виду изнашивания осуществлялось путём измерения поверхностной микротвердости на 
сканирующем нанотвердомере путём анализа толщины нанесенных индентором царапин. 

Результаты исследования. По результатам микроидентирования наибольшей  микротвердостью 
характеризуется экспериментальные сплавы 2.22 (связка 5,65%Со-1,8%Мо-+0,6%1Т1) и 2.23 (связка 
5,1%С0+2,7%Мо-+0,61%Т1). Для этих материалов средняя ширина царапины при различных усилиях была 
минимальна. В ходе трибологических испытаний наилучшие характеристики были зафиксированы для 
нержавеющей стали в сочетании со сплавом 2.22 и для пары трения «титановый сплав ВТЗ-1 — твёрдый сплав 
2.23». Процесс трения для этого сочетания материалов характеризуется невысокими коэффициентами трения с 
низким уровнем флуктуаций, минимальным износом образцов и изменениями начального микрорельефа их 
поверхностей. 

Обсуждение и заключения. В результате исследований установлены оптимальные с точки зрения 
трибологического взаимодействия пары трения: «титановый сплав ВТЗ-! — твёрдый сплав 2.23» и 
«нержавеющая сталь 12Х18Н9Т — твёрдый сплав 2.22». Процесс фрикционного взаимодействия для данных 
сочетаний материалов характеризуется минимальным объёмным износом, что будет способствовать 


повышению износостойкости инструмента на участках упругого контакта на передней и задней поверхностях. 
Ключевые слова: твёрдые сплавы, износостойкость, нержавеющая сталь, титановый сплав. 
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Введение. Различные материалы на основе карбида вольфрама, в том числе и твёрдые сплавы (ТС), 
получили широкое распространение во многих отраслях современного производства благодаря ряду 
преимуществ в области физико-механических свойств [1-3]. Наибольшая доля (порядка 65 %) применяемых в 
металлообработке инструментальных режущих материалов (ИРМ) также принадлежит ТС, обеспечивающим 
высокие скорости резания при обработке различных конструкционных материалов [1, 4]. В настоящее время 
развивается множество направлений, связанных с улучшением эксплуатационных характеристик ТС 
различными методами [5—7]. Одним из них является создание новых составов связок для карбидных фаз [8—10]. 
На кафедре «Металлорежущие станки и инструмент» Донского государственного технического университета 
разработаны экспериментальные твердые сплавы (ЭТС) на базе однокарбидного сплава ВК8 с различными 
типами модифицированных связок'. Несмотря на то, что кобальт благодаря своим свойствам является наиболее 
распространённой связкой для \!С, применение этого метала экономически невыгодно ввиду его дефицитности 
и высокой цены. Экспериментальные сплавы характеризуются большими значениями тепловой энтропии 
модифицированной кобальтовой связки (таблица 1) и, следовательно, меньшими термоэдс по отношению к 


карбиду вольфрама, что повышает электрохимическую устойчивость этих материалов [10-12]. Среди 


'Рыжкин А. А., Месхи Б. Ч., Боков А. И. и др. Твердый сплав на основе карбида вольфрама (варианты): патент 2531332: Рос. Федерация: 
МПК С22С29/08 № 2012128284/02: заявл. 04.07.2012, опубл. 20.10.2014. Бюл. №29. бс. 
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обрабатываемых с помощью твердосплавного инструмента материалов можно выделить титаново- 
аллюминиевые сплавы и аустенитные нержавеющие стали. За счёт ряда характеристик эти материалы 
применяются для изготовления деталей в наиболее ответственных областях машиностроения, включая 
аэрокосмическое, атомное, пищевое и медицинское производства. Обработка резанием титановых сплавов и 
нержавеющих сталей затруднена в связи с их низкой теплопроводностью, высокими усилиями при резании, а 
также неудовлетворительными триботехническими характеристиками [13]. В этой связи исследование 
особенностей фрикционного взаимодействия вновь разрабатываемых ИРМ с упомянутыми материалами 
является весьма важной актуальной задачей. Учитывая современные экологические и экономические 
требования к организации производства, обработка этих конструкционных материалов происходит 
преимущественно в режиме минимальной подачи (МОГ, — шипит даайщу Габисайоп) или без использования 
смазывающе-охлаждающих технологических сред (СОТС) [14, 15]. Тогда участки контакта ИРМ и материала 
обрабатываемой заготовки на передней и задней поверхностях режущего инструмента можно рассматривать 
как трибосистемы, функционирующие в режиме трения без СОТС или в режиме граничного трения. Это не 
исключает значительную долю металлического контакта. В этом случае на эксплуатационные параметры ИРМ, 
как элемента пары трения, будут оказывать значительное влияние его триботехнические характеристики в 
режиме сухого трения. 

Настоящее исследование посвящено изучению триботехнических характеристик экспериментальных 
твёрдых сплавов при сухом трении по труднообрабатываемым резанием материалам, а также определению 
оптимального с точки зрения фрикционного взаимодействия сочетания «ЭТС — конструкционный материал». 
Данная работа является частью комплекса исследований физико-механических, трибологических и режущих 
свойств экспериментальных ТС с модифицированной кобальтовой связкой. 

Материалы и методы. В качестве объектов исследования поверхностных механических характеристик 
выбраны нижеследующие составы ТС (таблица 1). Трибологические испытания проводились на инденторах 
квадратного сечения (а=5 мм, Ва=0,1-0,12 мкм) из трёх наиболее перспективных составов (2.21, 2.22 и 2.23) на 
трибометре Т-11(Польша), реализующего схему трения «палец-диск». В ходе проведения опытов 
фиксировалось изменение силы трения (Ё, Н) в зависимости от пути трения (Г, м). Каждый эксперимент 
повторялся 3—5 раз, результаты опытов обрабатывались при помощи методов дисперсионного анализа [16]. В 
качестве материала вращающегося диска (контртела) были выбраны титановый сплав ВТЗ-1 и нержавеющая 


сталь 12Х18Н9Т; шероховатость поверхности этих образцов находилась в пределах К. 0,12-0,15 мкм. 

































































Таблица 1 
Химический состав и термодинамические свойства ЭТС 
Энтропия, Дж/моль‘град 
Марка 
твёрдого Состав 
сплава ба ИС В связ. 
1 2 3 4 
92,63% \С-7,37%[1,52%Со-+ 
2.19 28,60 
+5,03%Ее+0,82%Си] ? 
2.20 92,38% \С-7,62% [3,6%Со-+ 29,31 
+3,2%Ее-+0,82%Си] 
2.21 92,45% \С+7,55% [5,3%Со+1,43%Ее+0,82%Си] 29,83 
2.22 91,95% \!С+8,05% [5,65%Со-+1,8%Мо-+0,6%ТИ 35,6 29,72 
2.23 91,59% \С+8,41% [5,1%Со-2,7%Мо-+0,61%Т 29,59 
2.24 90,62% \№\С+9,38% 29,22 
[3,34%Со-+5,44%Мо-+0,6%ТИ 
ВК8 92% УС - [7,5-8]%Со, Ее<0,3% 28,50 
(Базовый) 








Исследования проводились при различных скоростях скольжения и температурах при постоянной 


нагрузке Р=20 Н. Масса образцов определялась на весах ЛВ 210-А. Шероховатость дорожек трения на 
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контртеле после проведения экспериментов измерялась на профилометре Абл5-РМ7 (Россия). Определение 
динамической микротвёрдости поверхностей ЭТС производилось на сканирующем нанотвердомере 
Мапо5САМ-01 (Россия) путём анализа толщины нанесённых с различным усилием царапин. Исследования 
изношенных поверхностей образцов ЭТС осуществлялись на инвертированном микроскопе 7Е1$$ Ах1о Уем. 
АТ. 

Результаты исследования. Наименыпая ширина царапин й, нанесённых при различных усилиях Р., 


принадлежит сплавам 2.22 и 2.23, связка в которых была модифицирована группой Мо-Т! (таблица 2). 




















Таблица 2 
Ширина царапины й при различных усилиях А; по результатам склерометрии ЭТС 
й, мкм 
Е.Н Марки твёрдых сплавов 

2.19 2.20 2.21 2.22 2.23 2.24 ВК8 
5 0,3 0,2 0,5 - - 1.5 1,1 
15 1,3 1,3 1,5 0,7 1,2 2,5 1,8 
25 2,4 1,8 24 1,2 1,8 3,3 29 
































Таким образом эти материалы характеризуются наибольшей твёрдостью поверхности на микроуровне, 
что предполагают лучшую устойчивость к абразивному виду изнашивания. Следует отметить, что в этих 
сплавах при усилиях Е,;<15 Н заметных изменений поверхностного микрорельефа не было обнаружено. 
Наименьшую микротвёрдость по результатам испытаний продемонстрировал состав 2.24. 

Для определения и сравнения износостойкости ТС производилось измерение потери массы индентора 


для каждой величины пути трения Г, а затем определялся объёмный износ образцов ДТ (рис.1, 2). 
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Рис. 1. Объёмный износ инденторов ЭТС при трении по титановому сплаву ВТЗ-1 при температурах: 
а— 25°С; 6 — 300°С 
При трении по титановому сплаву при различных температурах наибольший объёмный износ ДИ 
наблюдался у сплава 2.22. Наилучшую износостойкость продемонстрировал состав 2.23 (рис. 1). При трении по 
нержавеющей стали наименьшие значения параметра ДУ были зафиксированы для состава 2.22 (рис. 2). 
Наибольший объёмный износ в этом случае также принадлежит экспериментальным составам ТС. При 
комнатной температуре самую высокую интенсивность изнашивания демонстрирует сплав 2.21, а при 


подогреве зоны трения максимальные значения износа фиксируются для состава 2.23. 
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Рис. 2. Объёмный износ инденторов ЭТС при трении по нержавеющей стали 12Х18Н9Т при температурах: 
а— 25°С; б — 300°С. 

Поверхности инденторов из сплавов, продемонстрировавших максимальный объемный износ, 
характеризуются либо преобладанием изношенных участков, либо изначальная поверхность материала 
сохраняется только в виде отдельных редких фрагментов. На рис.3 приведены снимки поверхностей 
инденторов из сплавов 2.23 и ВК8 после трения по стали 12Х18Н9Т при температуре 300°С для пути трения 
[=600 м. 





200 мкм 
а) 6) 


Рис. 3. Сравнение поверхностей инденторов сплавов 2.22 (а) и ВК8(б) 
после трения по нержавеющей стали 12Х18НЭТ (7=300°С, у=0,3 м/с): 
1 — изношенные участки поверхности; 2 — фрагменты исходной поверхности 


Поверхность более износостойкого сплава 2.22 имеет пятнистую структуру с преобладанием 
начальных неизношенных участков (рис. 3 а). Поверхность базового сплава ВК8 характеризуется большим 
масштабом разрушений, редкие фрагменты исходной поверхности расположены дискретно, их суммарная 
площадь при этом значительно меньше (рис. 3 6). 

Оценка фрикционного взаимодействия пар трения производилась путём сравнения средних значений 
коэффициента трения с и его флуктуаций без учёта стадии приработки 1 (рис. 4 а, 6). Сравнение колебаний 
коэффициента трения осуществлялось по среднеквадратическому отклонению ох» этого параметра от /‹» на 


стадии установившегося трения П. 
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Рис. 4. Зависимость коэффициентов трения Гот пути Г, при трении по сплаву ВТЗ-1 (Т=25°С, у=0,3 м/с): 
а — 2.23; в — ВКВ: [ — стадия приработки; П — стадия стабильного трения 





В случае трения по титановому сплаву при различных температурах средние коэффициенты трения и 
их среднеквадратичные отклонения у всех ЭТС были выше, чем для базовой марки ВК8. Однако, наименьшие 
значения этих параметров среди экспериментальных составов принадлежат сплаву 2.23 (рис. 4 а). 

При трении по стали 12Х18Н9Т как при комнатной температуре, так и с подогревом наибольшие 
значения параметров /.› и 0х, принадлежат базовому сплаву ВК8 (при 25°С: /,=0,72 и о. „=0,048; при 300°С: 
Лр=0,68 и 0‹.=0,032). Минимальные коэффициенты трения и его флуктуации в этой серии опытов были 
зафиксированы для состава 2.22. При трении без подогрева для этого материала значения оцениваемых 
параметров составили /.,=0,44 и о„=0,025, а при 300°С увеличились до /=0,57 и охр= 0,029. 

Произведено измерение и сравнение средних шероховатостей дорожек трения на контртелах для 
максимального значения пути трения. При трении по стали 12Х18Н9Т наименьшие значения этого параметра 
принадлежат сплаву 2.22. При температуре 25°С шероховатость составила Ка>›=4,12 мкм, при 300°С — 
Ка››2=5,12 мкм. Наибольшие значения шероховатости были зафиксированы при трении базового сплава ВК8. 
При комнатной температуре значение этого параметра составило Кавкз=5,07 мкм, с подогревом зоны трения — 
Кавкз=5,95 мкм. Высокие значения шероховатости свидетельствуют о том, что фрикционное взаимодействие 
этого материала с нержавеющей сталью сопровождалось более масштабными разрушениями в процессе 
образования и разрушением адгезионных и когезионных швов на поверхности материала контртела. 

В серии опытов с титановым сплавом ВТЗ-1 наилучшие показатели по шероховатости контртел были 
зафиксированы после трения экспериментального состава 2.23 (Ва›_з=3,35 мкм при 25°С и Ва›>з=4,54 мкм при 
300°С). Наибольшие поверхностные разрушения получили образцы после фрикционного взаимодействия со 
сплавом 2.22. Шероховатость поверхности при трении без подогрева составила Ва›2›=6,88 мкм, с подогревом 
— Ва) ›>2=8,07 мкм. 

Обсуждение и заключения. В результате исследования трибологических характеристик 
экспериментальных твёрдых сплавов при трении по труднообрабатываемым материалам установлены 
наилучшие сочетания с точки зрения фрикционного взаимодействия пары материалов. Для нержавеющей 
стали 12Х18Н9Т лучшие трибопоказатели были зафиксированы при трении в сочетании со сплавом 2.22, для 
титанового сплава — с составом 2.23. Процесс трения для данных сочетаний материалов при различных 
температурах характеризуется минимальным объёмным износом, невысокими коэффициентами трения и 
меньшим масштабом разрушений поверхностей как твёрдосплавных инденторов, так и контртел из 
конструкционных материалов. 
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Введение. Рассмотрены значение машинного обучения в условиях цифровой трансформации промышленности, 
методы реализации глубинного обучения для обеспечения производительности доверительных систем 
управления. Определена необходимость использования для глубинного машинного обучения сверточных 
искусственных нейронных сетей. Кратко рассмотрены различные технологии и архитектуры реализации 
искусственных нейронных сетей, проведен сравнительный анализ их производительности. Целью работы 
является исследование необходимости разработки новых подходов к архитектуре вычислительных машин для 
решения задач глубинного машинного обучения при реализации доверительной системы управления. 
Материалы и методы. В условиях цифровой трансформации использование искусственного интеллекта 
выходит на новый уровень. В основе технической реализации искусственных нейронных систем с глубинным 
машинным обучением лежит использование одной из трех базовых технологий: высокопроизводительных 
вычислений (НРС) с параллельной обработкой данных, нейроморфных вычислений (МС) и квантовых 
вычислений (ОС). 

Результаты исследования. Проведен анализ моделей реализации глубинного машинного обучения, базовых 
технологий и архитектуры вычислительных машин, а также требований по обеспечению доверия к системам 
управления, использующих глубинное машинное обучение. Выявлена проблема дефицита вычислительных 
мощностей для решения таких задач. Ни одна из существующих технологий не позволяет решать полный 
комплекс задач обучения и импеданса. Современный уровень технологий не обеспечивает информационной 
безопасности и надежности работы нейронных сетей. Практическая реализация доверенных систем управления 
с глубинным машинным обучением на базе имеющихся технологий для существенной части задач не 
обеспечивает достаточной производительности. 

Обсуждение и заключения. Проведенное исследование позволило выявить проблему дефицита 
вычислительных мощностей для решения задач глубинного машинного обучения. На основе анализа 
требований к доверенным системам управления определены объективные сложности их реализации на базе 
существующих технологий и установлена необходимость разработки новых подходов к архитектуре 
вычислительных машин. 


Ключевые слова: глубинное машинное обучение, процессор, доверенная система, информационная 
безопасность, вычислительная машина, искусственный интеллект. 
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Введение. В последние 9 лет в мире, и в России в т.ч., все более масштабно совершается 4-я 
промышленная революция, одной из ключевых составляющих которой является цифровая трансформация, 
затрагивающая все аспекты экономической жизни — от крупного промышленного производства до сферы 
обслуживания, науки, образования и домашних хозяйств. В этих условиях применение машины в задачах, с 
которыми человек не справляется или справляется хуже машины, существенно расширяется. Если раньше речь 
шла о механизации ручного труда и автоматизации производства, то теперь человек может быть заменен 
машиной при решении задач обработки данных, анализа, прогнозирования и управления различными 
системами: оборудованием, технологическими процессами, промышленными предприятиями, торговыми 
сетями и др. 

Практической основой замены человека машиной в отдельных областях интеллектуальной 
деятельности является использование искусственного интеллекта. Создание «сильного» искусственного 
интеллекта — задача будущего, связанная с необходимостью создания и развития новых технологий, а также 
решением значимых этических проблем. В настоящее время доступны для использования системы машинного 
обучения различной сложности, представляющие собой ступень к созданию «сильного» искусственного 
интеллекта. 
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Мировой рынок систем машинного обучения активно расширяется. В 2020 г. его объем составил 
11,33 млрд $, в 2021г. он вырос до 15,50 млрд $, ак 2028г. достигнет уровня 152,24 млрд $, показывая 
среднегодовой рост на 38,6 %'. 

Сфера применения систем машинного обучения очень велика и включает в себя маркетинг и торговлю, 
банковскую деятельность, промышленное производство, медицину и др. Системы машинного обучения 
наиболее востребованы в следующих областях промышленности: 

—щ в робототехнике для интеллектуализации промышленных и сервисных роботов, в т.ч. 
коллаборативных; 

— вавтоматизированных системах управления технологическими процессами и предприятиями; 

— в системах управления производственными процессами; 

— в системах управления поставками и взаимоотношениями с клиентами; 

— в исполнительных системах производства; 

— в системах производственной аналитики технологического оборудования; 

— всистемах бизнес-аналитики и др. 

Сложные технологические системы, например реализуемые в конструкциях металлорежущих станков с 
числовым программным управлением в режиме реального времени, в настоящее время не могут быть 
оснащены системами управления, пригодными для машинного обучения. Для этого требуются вычислительные 
мощности, способные выполнять существенные по объему вычисления в течении десятков микросекунд. Такая 
задача с помощью современных технических средств не может быть решена. Поэтому в большинстве случаев 
процесс машинного обучения системы управления осуществляется на отделенных от нее вычислительных 
машинах без ограничений времени, а затем результаты обучения передаются в систему управления в виде 
рекомендаций, указаний по режимам работы, смене инструмента, интервалам поверки ит. д. 

Методы машинного обучения условно соответствуют видам умозаключений, лежащих в их основе: 
индукции, дедукции и традукции. Метод контролируемого обучения по прецедентам, когда в машину 
загружаются большие объемы данных, предварительно маркированных оператором-человеком, соответствует 
индукции. Метод обучения без учителя, когда машина должна сама найти закономерности в данных, выявить 
шаблоны, упорядочить и структурировать данные, соответствует индукции и традукции. Дедукции и традукции 
соответствует экспертный метод, основанный на использовании для обработки данных заданных 
закономерностей и шаблонов. Традукция реализуется в основном через использование трансферного обучения, 
основанного на применении к данной задаче знаний, полученных при решении другой задачи. 

Для машинного обучения используются различные технологии и математические модели. Наибольший 
потенциал развития имеет модель искусственных нейронных сетей (ИНС), построенная по аналогии с 
биологическими нейронными сетями, т. е. сетями нервных клеток живого организма. ИНС представляет собой 
систему соединённых между собой и взаимодействующих искусственных нейронов, реализуемых в виде 
процессоров, процессорных элементов в виде ускорителей или сопроцессоров под управлением центрального 
процессора. Нейроны ИНС располагаются по уровням (слоям). Первый уровень соответствует получению, 
обработке входных данных и передаче их на следующий уровень. Промежуточные уровни — скрытые, их 
задача обрабатывать поступающие данные и передавать на последний (выходной) уровень. В нейронной сети 
может быть несколько скрытых уровней, перемежающихся с уровнями, на которых выполняются логические, 
математические и другие преобразования. От уровня к уровню данные обрабатываются, на каждом 
последующем уровне выявляются взаимосвязи предыдущего. Такая многоуровневая ИНС имеет большие 
возможности и может быть использована для реализации глубинного машинного обучения [1-3]. 

Глубинное машинное обучение — метод проектирования ИНС с помощью многослойных фильтров для 
извлечения и моделирования признаков из набора входных данных”. Такое обучение может быть 
контролируемым и неконтролируемым. Также возможно использование глубинного машинного обучения для 
экспертных систем. 

В основе технической реализации искусственных нейронных систем с глубинным машинным 
обучением лежит использование одной из трех базовых технологий: высокопроизводительных вычислений с 
параллельной обработкой данных, нейроморфных и квантовых вычислений°. 


МасЫпе Геагоше Маше 2021-2028 // Нага\маге & ЗоЁ\уаге ГГ Зегу1сез Маше Везеагсь Верог. — 2021. — Р. 160. — ОВС: 
Брз: ууу. Еогллие6 $ 1е551151265.сошЛи стары с /тасЬте-еагите-тагКе!-102226 (дата обращения: 06.11.2021). 
? Глек, П. Глубокое обучение (еер Тлеагпт?): краткий туториал // пеигоуело : [сайт]. ОВ: В рз://пеиго уе Ло/ли/озпоуу- 
Чавазслепсе/>абокое-обисвеше-4еер-Шеагиш2-Кта-виюона| (дата обращения: 06.11.2021). 
3 Как сократить издержки при использовании ИИ / НиасУатщагаСогроганоп : [сайт]. ОВ: Вирз://Виасв!.спемуз.пи/агс1е5/2021-06- 
14 КаК зокгайе 174ег7 1 ри _15ро!хоуапи (дата обращения: 07.11.2021). 
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Материалы и методы 

Высокопроизводительные вычисления. Высокопроизводительные вычисления с параллельной 
обработкой данных реализуются посредством гибридных вычислительных систем, т.е. систем с гетерогенной 
аппаратной вычислительной структурой, включающей в себя центральный процессор (СРО) и дополнительный 
вычислительный модуль в виде ускорителя или сопроцессора. В зависимости от процессоров, используемых для 
параллельной обработки данных, гибридные вычислительные машины имеют одну из четырех архитектур: 

1. Архитектура на базе графических процессоров (СОРИ). Наиболее распространенные решения 
представляют собой графические ускорители, расширяющие вычислительные возможности центрального 
процессора компьютерной системы. Последние достижения в данной области — графические ускорители МУШГА 
Теа \У100, обеспечивающие для задач глубинного машинного обучения производительность 120 ТЕГОРВ, т. е. 
1,2х1014 операций с плавающей запятой в секунду“. Это в 500-1000 раз выше производительности обычного 
персонального компьютера (ПК). Необходимо также учитывать, что указанная производительность обеспечивается 
при решении задач, требующих значительных вычислительных мощностей и существенных затрат времени, но не 
при работе в режиме реального времени. Архитектура на базе графических процессоров в настоящее время наиболее 
доступна. В частности, для реализации системы с ограниченной вычислительной мощностью достаточно иметь на 
ПК видеокарту с графическим процессором МУ, реализующую программно-аппаратную архитектуру 
параллельных вычислений СОПЛА. Наряду с СОРА к числу технологий СОРОРО, использующих графический 
процессор видеокарты для компьютерной графики в целях производства математических вычислений, относится 
технология АМР ЕиеЗиеат (для графических процессоров ускорителей АТП. Мировой рынок графических 
процессоров в настоящее время составляет около 26 млрд $, темп роста рынка — до 34 % в год°. 

2. Архитектура на базе программируемых пользователем вентильных матриц (ЕРСА) — 
полупроводниковых устройств, которые могут перепрограммироваться и менять топологию соединений в процессе 
использования. Номинальная производительность данных устройств сравнительно невысока — около 20 ТЕГОР$, 
однако эффективность использования вычислительных мощностей наибольшая среди всех рассмотренных 
архитектур. Она в 6-7 раз выше, чем у графических ускорителей. Высокая эффективность ЕРСА обусловлена 
гибкостью и скоростью подстройки под решаемые вычислительные задачи. Мировой рынок ЕРСА по данным 
компании Отап У1е\ Кезеагсв, в 2020 г. составил 9,85 млрд $, ожидаемые темпы роста рынка на период до 2027 г. 
— 97 %в годе. 

3. Архитектура на базе интегральных схем специального назначения (АЗС). Благодаря узкой 
специализации решаемых вычислительных задач они могут быть существенно проще, дешевле и компактнее. 
Производительность АЗС может достигать 1000 ТЕГОР$, однако эффективность использования вычислительных 
мощностей, например количество распознанных изображений, примерно в 2 раза ниже, чем в графических 
ускорителях. Темпы роста рынка АЗС существенно ниже, чем рынка графических процессоров. В 2020г. по 
данным компании С]оба| шдизгу Апа!у$ мировой рынок АЗС составил 17,3 млрд $, ожидаемый среднегодовой 
рост рынка до 2027 г. — 7,7 %'. 

4. Архитектура на базе однокристальных ускорителей (50С). «Система на кристалле» 50С — 
полнофункциональное электронное устройство, имеющее материнскую плату, процессор и другие необходимые для 
работы компоненты, размещённое на одной интегральной схеме. 5оС распространены в мобильных компьютерах 
(смартфонах), одноплатных компьютерах и других встроенных системах. При этом $0оС имеют существенный 
потенциал использования в составе гибридных вычислительных машин. Кроме того, возможны решения в виде 
однокристальной сборки З0С с элементами ЕРСА (архитектура Уегза| для адаптивных вычислений компании 
Хих°). Рынок $0С в настоящее время очень велик и за 2020г. составляет 79,7 млрд $. На период до 2027 г. 
прогнозируется рост рынка на 4,4 % в год, к 2027 г. объем рынка достигнет 107,4 млрд $°. 

На рис. 1, 2 приведены данные по фактическому и прогнозируемому росту мирового рынка чипов для 
глубинного машинного обучения, подготовленные консалтинговой компанией Оша °°. 


4 Как сократить издержки при использовании ИИ Указ. соч. 
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Рис. 1. Динамика роста количества чипов для глубинного обучения по годам 
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Рис. 2. Динамика изменения рынка чипов для глубинного обучения по годам"! 


Анализ показывает, что на фоне общего роста рынка наибольшие перспективы имеют АЗС, значимые 
позиции сохранят также СРО и 50С. Ускорители СРИ на достаточно длительную перспективу не имеют 
адекватной замены в решении сложных задач, в т. ч. в процессе обучения, а ЗоС незаменимы для мобильных 
реализаций систем глубинного машинного обучения, а также для параллельных вычислений с целью разгрузки 
центрального процессора. Отказ от ЕРОА в системах глубинного машинного обучения представляется 
необязательным, хотя доля таких процессоров, вероятно, будет сравнительно небольшой. 

Одним из наиболее перспективных направлений развития гибридных вычислительных машин является 
использование тензорных и других специализированных сопроцессоров, например процессоров машинного 
зрения. Такие сопроцессоры могут строиться на базе наиболее распространенных и производительных АЗС, а 
также ЕРСА или СРО. Отличие сопроцессоров от ускорителей заключается в степени интеграции с 
центральным процессором. Центральный процессор посредством специальной области памяти транслирует 
управляющие инструкции на ускоритель. Сопроцессор отслеживает поток инструкций машинного кода, 
поступающий из оперативной памяти в центральный процессор, и перехватывает инструкции, 
соответствующие его функциональному назначению, например задачи тензорных преобразований, 
распознавания образов и т.д. Для решения масштабных задач, требующих длительных распределенных 
вычислений, целесообразно применение ускорителя, для частого и многократного выполнения несложных 
вычислительных задач — применение сопроцессора, при котором центральный процессор не загружается и не 
замедляет обработку данных. 

Эксплуатационные свойства вычислительных машин существенно зависят от используемой 
архитектуры. Для решения задач глубинного машинного обучения, куда затрачиваются до 80% 
вычислительных мощностей, наилучшим образом подходят машины на базе графических ускорителей. Они 
имеют высокую производительность при решении сложных задач, требующих существенных затрат времени, 
обладают высокой гибкостью и максимальной точностью вычислений, но низкой относительной 
производительностью. Например, для процессоров МУ1Ча она составляет 1,3-1,8 СОРУ/\М/. Вычислительные 


И ТозЪЬ, А. Оеер Геагите СЫ рзев Вером — 2020 // От@а Маше. — 2020. — ОВГ: Вирз://ота!алесВ.{поггла.сога/ргоисЁ5/еер-еатте- 
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машины на базе АЗГС обладают максимальной абсолютной и относительной производительностью. Например, 
для процессоров пеВМ относительная производительность составляет 254 СОРЗ/Л\ [4]. Однако такие 
машины имеют низкую гибкость и ограниченную точность, поэтому их целесообразно использовать при 
решении типовых, например матричных или тензорных преобразований, повторяющихся или многопоточных 
задач в режиме реального времени. 

Вычислительные машины на базе ЕРОА обладают высокими параметрами гибкости, точности, 
абсолютной и относительной производительности. Например, для процессоров ТезтаТХ1 относительная 
производительность составляет 70 СОРУ/\ Однако такие машины имеют сравнительно высокую стоимость, 
поэтому их целесообразно использовать для научных целей, когда требуется лишь несколько вычислительных 
машин заданной конфигурации, а также для разработки архитектуры массово выпускаемых процессоров АЗС 
и 50С. 

Несмотря на совершенствование архитектуры вычислительных машин потенциал роста мощностей для 
высокопроизводительных вычислений в скором времени будет исчерпан. Число транзисторов на кристалле за 
последние 5 лет увеличилось примерно в 12 раз [5], а объем вычислений в процессе машинного обучения — в 
150 тыс. раз". 

Нейроморфные вычисления. Возможным способом устранения нехватки вычислительных 
мощностей для искусственных нейронных систем с глубинным машинным обучением является использование 
нейроморфных вычислений и соответствующих чипов. Нейроморфный чип — это процессор, работа которого 
основана на принципах человеческого мозга. Такой чип моделируют работу нейронов и их отростков — 
аксонов и дендритов, отвечающих за передачу и восприятие данных. Связи между нейронами образуются за 
счет синапсов — специальных контактов, по которым транслируются электрические сигналы. 

К числу наиболее известных разработок в данной области относятся нейроморфные процессоры 
ТгаеМойй компании ВМ и процессоры Го! компании ш@]. В них используется асинхронная кластерная 
архитектура и модель сверточной нейронной сети — однонаправленная многослойная сеть с чередованием 
свёрточных и субдискретизирующих слоёв. Процессор ТгаеМогй выполнен на базе технологий 28 нм, Го — 
на базе 14 нм [6]. 

Многопроцессорная система ТгиеМонн №16е-4, содержащая 100 млн нейронов и предназначенная для 
работы с сетями для глубинного машинного обучения, была представлена компанией 1ВМ в 2018 т. [7]. Каждый 
чип содержит 1 млн цифровых нейронов и 256 млн синапсов, заключенных в 4096 синапсных ядрах; 
энергопотребление каждого чипа — 70 мВт. 

Процессор Гой\, представленный в 2017 г., содержит 131 тыс. искусственных нейронов и 131 млн 
синапсов. В 2019-2020 гг. Пи представила два продукта на основе Гот — Робой\ВеасВ и Рово$рииз$. 
Вычислительная система РоройчВеасв, включающая 64 процессора Гот, суммарно располагает 8,32 млн 
нейронов и 8,32 млрд синапсов. Вычислительная система Рово 1$ рипя$ включает 768 процессоров Го, имеет 
100 млн нейронов и 100 млрд синапсов". 

В России работы по созданию нейроморфных процессоров ведутся уже несколько лет. В 2020 г. 
компания «Мотив Нейроморфные Технологии» создала нейрочип «Алтай» [8]. Технологическая норма 
процессора — 28 нм, энергопотребление — около 0,5 Вт, площадь кристалла — 64 мм” (для сравнения: 
ТгиеМогй — 430 мм”, Го! — 60 мм). В нем 131 тыс. нейронов, между ними — 67 млн связей. 

Для оценки качества нейроморфных процессоров используют: 

1. Абсолютный показатель производительности. Это количество миллиардов выполняемых 
синоптических операций в секунду (СЗОР5). 

2. Показатель энергоэффективности. Это количество пикоджоулей энергии, затраченной на 
выполнение одной осинаптической операции (р//ЗОР). 

Процессор ТгаеМойй имеет производительность 58 9$ОР$ и энергоэффективность 26 р//ЗОР\. 
Аналогичную энергоэффективность (23,7 р/ЗОР [9]) имеет процессор Гойи. 

Единственным конкурентом нефроморфных процессоров в реализации нейронных сетей с глубинным 
машинным обучением на среднесрочную перспективу (8—12 лет) являются гибридные вычислительные машины 
с сопроцессорами АЗС. Такие процессоры имеют меньшую, но сопоставимую производительность 
выполнения синаптических операций и более высокую энергоэффективность. В частности, процессор АЗС, 





йе ТВоштрзоп, М.С., Отеепе\ууа]а, К., Гее, К., Мапзо, С.Е. ТВе Сотриайопа| ГлтИз оЁОеер Геагипе// агХлу ргерии агХау:2007.05558. — 2020. 
(дата обращения: 12.11.2021) 

3 пи& Зса!ез МеиготогрЫс Везеагсв бузет 10 100 МИНоп Мешгопз// пие|. — 2020. — ОВГ: Вирз://пе\узгоота. не. соталлезуз/либе]-зса[ез- 
пеиготогрЫс-гезеагсв-зузет-100-пИПоп-пеигоп5/ (дата обращения: 12.11.2021). 

1 Нейроморфный процессор «Алтай» // тоНупети : [сайт]. ОВГ: Врз:/лпоНупги/пеигосыр-аЦа!/ (дата обращения: 10.11.2021). 
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описанный в работе [10], обеспечивает производительность выполнения синаптических операций 8,7 СЗОР$ и 
энергоэффективность 15,2 р/ЗОР. 

Мировой рынок нейроморфных чипов сформировался сравнительно недавно и поэтому невелик. В 
2020 г. его объем составил всего 22,5 млн $. При этом темпы роста рынка очень высоки. К 2026 г. рынок 
вырастет до 333,6 млн $, что соответствует среднегодовому росту на 47,4 %”. 

Иногда под нейроморфными чипами понимают все виды процессоров, внешне воспроизводящих 
работу нейронов независимо от внутренней структуры технического устройства, которое может не 
соответствовать характеру взаимодействия нейронов. Такие процессоры, используемые для построения 
искусственных нейронных сетей, правильно называть нейронными. К нейронным процессорам, наряду с 
нейроморфными чипами, также относятся процессоры (чипы) с тензорными и другими специализированными 
сопроцессорами машинного зрения, распознавания речи и др. Мировой рынок нейронных процессоров в 
настоящее время составляет 2,3 млрд $, ак 2027 г. вырастет до 10,4 млрд $, т. о. средний рост составит 24,2 % в 
годе. 

Квантовые вычисления. В долгосрочной перспективе средством устранения дефицита 
вычислительных мощностей может стать развитие квантовых вычислений. Квантовые вычисления — решение 
задач с помощью манипуляции квантовыми объектами: атомами, молекулами, фотонами, электронами и 
специально созданными макроструктурами. Манипуляции квантовыми объектами дают возможность 
использования: 

—щ квантовой суперпозиции, которая проявляется в способности квантовых систем одновременно 
находиться во всех возможных состояниях; 

—щ квантовой запутанности, которая проявляется в сильной взаимосвязи между параметрами 
специальным образом приготовленных квантовых систем. 

Устройства для квантовых вычислений принято разделять на два больших класса [11]: универсальные 
квантовые компьютеры и квантовые симуляторы. Первые подобно центральным процессорам могут решать 
любую алгоритмическую задачу, квантовые симуляторы — это аналоговые компьютеры для решения 
узкоспециализированных задач. 

Технологии создания универсальных квантовых компьютеров в настоящее время находятся на стадии 
формирования. Создаваемые вычислительные машины демонстрируют «квантовое превосходство» в решении 
отдельных задач, но пока не могут быть использованы для формирования искусственных нейронных сетей с 
глубинным машинным обучением. Наибольшую активность в создании квантового компьютера демонстрируют 
компании: 

1. Соозе. В 2018 г. построен 72-кубитный квантовый процессор ВизЧесопе, в 2019 г. — более точный 
53-кубитный квантовый процессор Зусатоге. 

2. Пие|. В 2018 г. построена 49-кубитная сверхпроводящая квантовая микросхема ТапеГаКе. 

3. ВМ. В 2017г. создан и протестирован 50-кубитный квантовый процессор, в 2019 г. — первый в 
мире коммерческий 20-кубитный квантовый компьютер 1ВМ О ЗузетОпе и др. 

Единственным адиабатическим квантовым компьютером, представленным на рынке, является О- 
\\!ауебузет$, выпускаемый в вариантах от 16 до 2000 кубитов, организованных в кластеры по 8 кубитов в 
каждом. 

Область квантовых симуляторов также активно развивается. Один из наиболее сложных симуляторов 
такого типа представлен совместной разработкой 2017г. Мэрилендского университета и Национального 
института стандартов и технологий (США). Этот 53-кубитный симулятор использует в качестве кубитов 
холодные ионы иттербия. Аналогичный по возможностям 51-кубитный квантовый симулятор на базе атомов 
рубидия был разработан группой ученых Гарвардского университета и Массачусетского технологического 
института. 

В России также реализуется ряд проектов, развивающих технологии квантовых вычислений. В 
частности, уже несколько лет ведётся разработка сверхпроводящего процессора учёными консорциума, куда 
входят МИСиС, ИФТТРАН, МГТУ им. Н. Э. Баумана, ВНИИА им. Духова и другие организации. На 
сегодняшний день консорциумом отлажена технология изготовления сверхпроводниковых двухкубитных схем, 


'° Мепготогрыс СЫр Маше — Ого\ж®, Тгепдз, СОУТО-19 Ипрасё апа Еогесазз (2021-2026) // Могаог ПиеШеепсе. — 2020. — ВТ: 
Брз://уумуу поогаогииеШоепсе.сот/а4изу-геройз/пеиготогрЫс-сЫр-пагКе! (дата обращения: 11.11.2021). 

° Мешготогрыс СЫрз — СОюба:! Маше Тгадесюгу & Апауйсз // Соба! шаизму Апайузз, ше. 2021. 118 р. ВЕ: 
Би рз:/Лууилу гезеагсвап4тагКе.сот/герог5/4805280/пеитотогрЫс-сЫрз-вюБа]-тагке-га]есюгу-ап4 (дата обращения: 10.11.2021). 
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экспериментально охарактеризованы и  продемонстрированы двухкубитные логические вентили, 
осуществляющие квантовое запутывание, необходимое для работы квантовых процессоров. Достоверность 
логических операций находится в пределах 85—95 %. 

В 2020 г. в квантовые вычисления в мире было вложено 675 млн $, что более чем в 3 раза превышает 
объем инвестиций 2019г. (211 млн $). В 2021г. объем инвестиций в квантовые вычисления превышает 
800 млн $ [12]. 

Доверительные системы управления. Одним из базовых требований к системам управления, в том 
числе к производственным, является обеспечение к ним необходимого уровня доверия. Согласно ГОСТ Р 
54583-2011 «Методы и средства обеспечения безопасности. Основы доверия к безопасности информационных 
технологий. Часть 3 — Анализ методов доверия» целью обеспечения доверия является создание уверенности в 
надёжном функционировании продукта в заданных условиях. Для обеспечения указанного доверия 
информационная система должна обладать следующими эксплуатационными свойствами [13]: 

—щ функциональной надёжностью, т.е. способностью выполнять свою функцию с заданной 
достоверностью, которая в свою очередь нормируется числом отказов, погрешностью и повторяемостью 
результатов вычислений; 

— информационной безопасностью, т.е. способностью обеспечивать заданный уровень 
конфиденциальности, доступности и целостности информации: хранимой, передаваемой, принимаемой и 
обрабатываемой в процессе работы системы. 

Предметом настоящей работы являются системы управления, которые относятся к информационным. 
Поэтому для них также справедливы указанные выше требования к эксплуатационным свойствам. Однако у 
систем управления имеется своя специфика в определении доверия. Доверенная система управления должна 
обладать: 

— способностью управлять, например роботом, станком, предприятием и т. д. по заданному числу 
параметров с заданной надежностью, которая регламентируется числом отказов, погрешностью и 
повторяемостью, и с заданной производительностью, которая в свою очередь регламентируется временем 
обработки данных и выполнения управляющих команд; 

— возможностью контроля элементов, структуры и процессов системы на аппаратном и программном 
уровне с целью обеспечения информационной безопасности. 

Если система управления имеет функцию глубинного машинного обучения, то выполнение первого из 
указанных требований накладывает жесткие ограничения на используемые средства технической реализации. 
Это должна быть оптимальная для формирования сверточной нейронной сети вычислительная машина с 
высокими параметрами производительности, точности и надежность вычислений. 

Если не рассматривать вариант использования квантовых процессоров, полнофункциональных 
реализаций которых пока нет, то первому требованию в полной мере не соответствуют ни гибридные 
вычислительные машины на базе всех рассмотренных архитектур, ни вычислительные машины на базе 
нейроморфных процессоров. Вычислительные машины на базе АЗС и нейроморфных процессоров не 
обеспечивают высокой точности и надежности, а гибридные вычислительные машины с ускорителями СОРИ или 
ЗоС неоптимальны для работы в режиме реального времени, в том числе импеданса. 

Определенный компромисс обеспечивается при использовании гибридных вычислительных машин на 
базе ЕРСА, однако такие машины имеют высокую стоимость при серийном производстве, существенно 
меньшую производительность, чем машины с АЗС, и существенно меньшие возможности сложных 
вычислений, чем машины с СРО. Другим вариантом достижения компромисса является одновременное 
использование центрального процессора с графическим ускорителем для решения сложных задач в процессе 
машинного обучения и тензорных или иных узкоспециализированных сопроцессоров на базе АЗС для 
обработки данных в режиме реального времени. 

Второе требование, хотя и является техническим по содержанию, на практике выступает в качестве 
экономического. Реализация системы управления технологическим оборудованием с глубинным машинным 
обучением возможна только посредством сверточной нейронной сети, контроль работы которой извне не 
представляется возможным. Информационная безопасность может быть обеспечена только при условии, что 
основная часть вычислительной машины будет создаваться отечественными производителями, которые прошли 
сертификацию в области информационной безопасности. 

В настоящее время в России основная часть систем управления построена на основе зарубежных 
микроэлектронных компонентов. Доля таких компонентов превышает 85 % [14]. Обеспечение информационной 
безопасности в случае использования в вычислительных машинах импортных компонентов не имеет 
однозначного решения и зависит от структуры создаваемой искусственной нейронной сети и порядка ее 
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использования. В частности, при использовании гибридных вычислительных машин информационная 
безопасность может быть существенно повышена за счет локализации передачи данных между центральным 
процессором и ускорителем или сопроцессором. 

Результаты исследования. Проведенный анализ моделей глубинного машинного обучения, базовых 
технологий и архитектуры вычислительных машин, а также требований по обеспечению доверия к системам 
управления, использующих глубинное машинное обучение, позволяет сделать следующие выводы: 

1. Имеется объективная проблема дефицита вычислительных мощностей для решения задач 
глубинного машинного обучения. Ни одна из существующих в настоящее время технологий не позволяет 
решать полный комплекс задач обучения и импеданса. 

2. Поскольку глубинное машинное обучение реализуется на базе модели сверточных нейронных сетей, 
их внешний контроль с целью обеспечения информационной безопасности и надежности работы не 
представляется возможным. Единственным вариантом является контроль разработчика, который также имеет 
ограниченные возможности. Это определяет необходимость производства необходимых для ИНС процессоров 
в России. 

3. Практическая реализация доверенных систем управления с глубинным машинным обучением на базе 
имеющихся технологий для существенной части задач в режиме реального времени не может быть обеспечена, 
для другой части задач такая реализация связана с существенным падением производительности. 

4. Повышение производительности доверенных систем управления может быть основано на 
совершенствовании архитектуры гибридных вычислительных машин, в том числе одновременном 
использовании процессоров разной архитектуры, оптимальных для решения соответствующих задач анализа и 
управления. 

Обсуждение и заключения. В настоящей работе проведен анализ и рассмотрены вопросы значения и 
реализации машинного обучения в условиях цифровой трансформации промышленности. Затронутая в работе 
научная проблема заключается в недостаточном развитии технического уровня современных вычислительных 
машин для обеспечения высокой производительности алгоритмов на основе глубинного машинного обучения. 
Обращает на себя внимание проблема обеспечения информационной безопасности, что является одной из 
предпосылок развития отечественных процессоров для ИНС. На основе анализа требований к доверенным 
системам управления определены объективные сложности их реализации на базе существующих технологий и 
необходимость разработки новых подходов к архитектуре вычислительных машин. 


Библиографический список 

1: От пеиго-ассеегабог$ оп ЕРСА$ ш соПабогайуе тобойс$ {азКз / 
А. ХаепзКу, Е. ЗетегзВ свету, А. МерКо [е# а1.] // ЗРТЕ Орйса! пзилитепе Зс1епсе, ТесБпо]оэу, апа Арр|Нсайопз П. 
— 2021. — \Уа/. 11876. — Ан. 1187600. — Р. 5. В рз://401.оге/10.1117/12.2600582 

2. ГеепзКи, А. А. Сопёго] оЁ СоПабогайуе Кобо{ Зузетаз апа Еех1Ые Ргодисйоп СеП$ оп Ше Ваз1$ о Оеер 
Геагише / А. А. Х@епзКи, М. М. Р1зтепзКота, У. У. Уогопи // Киз ап Епешеегте Везеагсв. — 2019. — У\Уо1. 39. 
— Р. 1065—1068. 6 рз://401.0г2/10.3103/5$1068798Х19120256 

3. Ачющтае4 у15иа| шзресйоп оЁ абмс ппазе изше 4еер Пеагите арргоасВ Юг 4еесё деесйоп / 
У. У. Уогопи, В. А. З1хуаКш, М. 7Вдапота [е1 а1.] // Ааютаеа У1зиа! шзресйоп ап Масбше У151юп ТУ. — 2021. 
— \Уо/. 11787. —Р. 117870. В рз://401.оге/10.1117/12.2592872 

4. Ме оу: РаёаЙо\ У151оп Ргосеззте Зузет-оп-а-Сыр / Ры1-Нипе Р|ват, О. Лейаса, С. Еагафей [е{ а1.] // 
шт: Ргос. ТЕЕЕ 55 Пиегпайопа! М14\ез{ бЗутрозциа оп Сисий5 апа Зузетз (МУ/ЗСА$). — 2012. — Р. 1044— 
1047. В@рз://401.оге/10.1109/М\У$СА$.2012.6292202 

5. Шуремов, Е. Л. Стоит ли увлекаться Большими Данными? / Е. Л. Шуремов // Учет. Анализ. Аудит. 
— 2020. — Т. 7, № 2. — С.17-29. Вирз://401.0т®/10.26794/2408-9303-2020-7-2-17-29 

6. То\уагаз агайс1а] сепега| пиеШоепсе \уиф вубма Тлапис сЫр агсЬесвте / Лиз Ре, Ге! Оепз, Зеп З0п5 
[её а|.] // Мабаге. — 2019. — Уо/|. 572. — Р. 106-111. Б&рз://401.0г2/10.1038/341586-019-1424-8 

7. Модха, Д. ТгаеМойй: от нуля к 64 миллионам нейронов / Д. Модха // Открытые системы. СУБД. — 
2019. — №3. — С. 8. 

8. ТгаеМомВ: 4ез1еп апа 00| Но\ оРа 65 ш\у 1 шИШоп пеигоп ргоэтатта е пеигозупарис сШр / 
Е. АКоруап, /]. Зауа4а, А. Саззау [ей а1.] // ТЕЕЕ ‘гапзасйоп$ оп сотршег-а14е4 4ез1еп оЁ ищеэтайе4 стсийз апа 
зузетв. — 2015. — \Уо/|. 34. — Р. 1537-1557. 6 рз://401.ог2/10.1109/ТСАО.2015.2474396 

9. Гот: А пеиготогрЫс тапусоге ргосеззог у оп-сЬар 1еагите / Мще Рау1ез, Магауап Згиуаза, Тзипэ- 
Нап Глп [её а1.] // ТЕЕЕ Мпсго. — 2018. — Уо]. 38. — Р. 82-99. № рз://401.0г2/10.1109/ММ.2018.112130359 


























Информатика, вычислительная техника и управление 


с^ 
л 


Бйр://уезиК-Чоп$а.га 


А4уапсей Епятеетия Везеагсй 2022. Т. 22, № 1. С. 57—66. 155М№ 2687—1653 





10. ЕРйслепе зупарзе тетогу з#тасеиге Юг гесопйзига Ме Ча! пеиготогрЫс Бага\уаге / 7. Кит, /. Коо, 
Т. Ки, .. Г. Кип // ЕгопЧег$ ш пеигозсепсе. — 2018. — У]. 12. — Р. 829. В рз://401.0г2/10.3389/пз.2018.00829 

11. Федоров, А. Квантовые вычисления: от науки к приложениям / А. Федоров // Открытые системы. 
СУБД. — 2019. — №3. — С. 14. 

12. \Ба{ Нарреп$ \Веп ТР Тиглз ю ‘\УБеп’ ш Оцапат Сотшрийпе? / /. Е. ВоЫег, М. Гапелопе, Е. Тао [ей 
а1.] // ВСС Плюца| ТгапзРЮогтайоп. — 2021. — Р. 20. 

13. Сабанов, А. Г. Доверенные системы как средство противодействия киберугрозам / А. Г. Сабанов // 
Защита информации. Инсайд. — 2015. — №3 (63). — С. 17-21. 

14. Каляев, И. А. Доверенные системы управления / И. А. Каляев, Э. В. Мельник // Мехатроника, 
автоматизация, управление. — 2021. — Т. 22, № 5. — С. 227-236. № рз://401.ог2/10.17587/тач.22.227-236 











Поступила в редакцию 27.12.2021 
Поступила после рецензирования 17.01.2022 
Принята к публикации 18.01.2022 


Об авторах: 

Зеленский Александр Александрович, директор Института цифровых интеллектуальных систем 
Московского государственного технологического университета «СТАНКИН» (127055, РФ, г. Москва, 
пер. Вадковский, 3 а), кандидат технических наук, доцент, Зсориз, Кезеагсвег, ОВ СТО, ИеепзКуаа@отайЙ.сот 





Абдуллин Тагир Хабибович, преподаватель кафедры «Промышленная электроника и 
интеллектуальные цифровые системы» Московского государственного технологического университета 
«СТАНКИН» (127055, РФ, г.Москва, пер. Вадковский, 3 а), ведущий инженер, бсориз, ОВС, 
еуегезбаитае(@уап4ех.ги 








Жданова Марина Михайловна, младший научный сотрудник ФГБОУ ВО «Московский 
государственный технологический университет «СТАНКИН» (127055, РФ, г. Москва, пер. Вадковский 
переулок, 3 а), Зсориз, Везеагспег, ОК СТО, тшр1зтепзКоуа(@тай.ги 





Воронин Вячеслав Владимирович, заместитель директора Центра когнитивных технологий и 
машинного зрения Московского государственного технологического университета «СТАНКИН» (127055, РФ, 
г. Москва, пер. Вадковский, За), кандидат технических наук, доцент, 5сориз, Кезеагспег, ОКСО, 





уогопш_$1@тай.га 





Грибков Андрей Армович, директор Аналитического центра Московского государственного 
технологического университета «СТАНКИНУ» (127055, РФ, г. Москва, пер. Вадковский, За), доктор технических 
наук, профессор, Зсориз, Везеагсвег, ОВ СТО, ап4агто(@уапаех.ги 





Заявленный вклад соавторов: 

А. А. Зеленский — формирование основной концепции, цели и задачи исследования, подготовка 
текста, формирование выводов; Т. Х. Абдуллин, М. М. Жданова — проведение исследования, анализ 
существующих подходов; В. В. Воронин, А. А. Грибков — анализ результатов исследований, доработка текста, 
корректировка выводов. 


Все авторы прочитали и одобрили окончательный вариант рукописи. 


А4уапсей Епятеетия Кезеагсй 2022. Т. 22, № 1. С. 67—75. 155М№ 2687—1653 





ИНФОРМАТИКА, ВЫЧИСЛИТЕЛЬНАЯ 
ТЕХНИКА И УПРАВЛЕНИЕ 
1УМЕОВМАТТОМ ТЕСНМОГОСУ, СОМРОТЕВ 
ЭСТЕМСЕ, АХО МАМАСЕМЕМТ 





т СНеск ог ирамез 
ООС 004.11.5 Опета[ агие 


№ рз://401.0г2/10.23947/2687-1653-2022-22-1-67-75 
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Вгои Копаше Ато5 © 2-3, Туап Зшйпоу © ‚ МабоиВ Мое Негтапп © [=]. 
Реор|ез' Ейепа$ р ОтуетзИу оЁ Виза (КООМ) (Мозсох,, Киззап Ееаеганоп) ” 
Ч ыгоцконатеатоз9@етай.сот 





Тпго4дисйоп. Согопауилз, а1з0 Кпо\ аз СОУТ-19, \уаз Ягз{ деесеа ш У/Гавап, СШпа, ш ОесетьЬег 2019. [ 1за 
атПу оРуп$е$ гапеште Нот Фе соттоп со] ю зеуеге асше гезриаюгу зупаготе (ЗАВ 5). ТВе зутрюил$ оЁзисй а 
уз аге зпиПаг 1ю фозе оРа со!4 ог зеазопа| аПеголез. ГЛКе ойфег гезриаюгу уга$ез, И 15 ташу тапзпийед тгоцэВ 
апбогпе 4гор|е{$ \Пеп соизШпя ог зпее7те. ТвегеЮге, е гесоспоп оР СОУТ-19 тедишез саге 1абогафогу 
апа]уз1$, ап Фе тедисНоп оР гесост оп гезоигсез 15 а та]ог сваПепзе. Оп 11 МагсВ, 2020, фе У\Уопа Неа 
Огоаттаноп (\УНО) 4ес1агеа СОУТО-19, саазе4 Бу ЗАВ $-СоУ-2, а рап4епис, аз Шеге Ва Бееп ап ехропепна| 
шсгеазе ш сазез уо|4\ме, ап4 4етапа юг пмепзуе Без апа ге!а{е4 зеасвигез Ва Раг ехсее4е4 ех15Нпз сарасИу. 
Тве ЙЯгз( ехаплр!ез ог 15 аге 1е геолопз оЁ НаЙу. ВтатИ гезлйеге@ фе Йгзё сазе ор ЗАВ$ -СоУ-2 оп 02/26/2020. 
Тгапзил1$1оп о# Фе ута$ ш 15 соипу зы Йе4 уегу диюеЮу Йот ппропе4 сазез © 1оса| апа, АпаПу, сотилипИу 
111551018, \УИВ фе Вта’Шап Ее4ега| соуегитепе аппоичпсте пайопа| соттипйу гапзпл15$10п оп 03/20/2020. Аз оЁ 
Магсь 23, ш Фе 5ае оЁ 5ао Рашо \И а роршайоп оЁ абоиё 12 шИПоп реоре, \Беге Бе [згаеШа АШег Етует 
Нозриа! 15 1осаёе4, 477 сазез оЁ {те @зеазе ап4 30 ге е4 4еа 5$ \еге геслуеге4, ап оп МагсН 27, феге \уеге атеа4у 
1223 сазез оГ СОУШ-19 уу 68 сопсотИапЕ 4еа $. То $10\и Ше зргеа оЁ Ше упа$ ш фе 5ае оЁ ЗАо Рау, 
апагапипез ап4 зос1а1 Из${апсте теазигез \еге шего4исе4. Опе оЁе тоНуаНопз$ Гог 1$ сраПепее 1$ {1е Ёас{ Фаф, ш 
Фе сощехЕ оРап ежепзуе Беа саге зузет \/ИВ Фе роз ]е шийайоп оЁ ЗАВ $-Со\У-2 {езИпо, И 15 по{ ргасйса! №0 
{е5{ еуегу сазе, ап4 {е5ё гези5$ сап ошу Бе изе4 ш 1езЯпе Фе {агоей заброршШайоп. ТВе заду обесйуе 15 0 БаЙ4 а 
то4е! Базе оп тасЬше |еагппе фа{ сап рге@с! фе деесйоп оГ ЗАВ 5$ -СоУ-2 Нот те1са| дайа. Еог 1$, уапои$ 
с1а551Исайоп по4е!5 оР тасЬше |еагите аге сотраге4, ап е Без( опе №ю рге41с( согопаупта$ез 1$ дейеглтеча. ТБе 
сотрагзоп 1$ Базе4 оп шт\у1Яча!$ т с]аз$ 1, 1.е., озе \уИВ а роз1йуе {езё. ТВегеге, й 15 гедите4 10 дейегите те 
тасЬше |еагипе то4е| у Бе Без гезропзе ап@ Е1 зсоге ог <1аз$ 1. 

Маептаб апа Мейоа$. Ап ореп-зоигсе дайа зе{ Нот Фе 15таеШа Абем Ешует Нозриа! ш Зао Рацо, Вга71, \’аз 
{аКеп а а Раз1з. Тре ЮПо\уше тасБше 1еагпте то4е]з ууеге изе4 Юг Фе зшау: КапдотРогез5 (КЕ), К-Меагез 
Метойбог (КММ), биррог Уесюг МасЬше (ЗУМ), Го215Яс Весотезз1оп (Г.В), Рес1з1оп Тгее (ОТ) апа АдаВооз (АВ), 
аз \е| аз Ше 10-ите сгоз5-уа|4айоп 1есЬмаце. Зоте тасЬше 1еагише регРоптапсе теазигез, засВ аз ассигасу, 
теса|, ап4 Е1 зсоге \уеге еуааае4. 

Кезий$. Ош оРа 1юа1 о 5,644 реоре 1е5$е4а 4игте Фе СОУТ-19 рапаепис, 5,086 реоре 1е5%е4 пегануе апа 558 
реор!е 1е$е4 розшуе. А Ше зате Ите, зиррог Рог тасЫте уесюгз зпо\е4 Ше Без гези $ ш 4еесИипе согопауйа$ 
УЕ а гесаП о 75 % ап4 ап Е1 зсоге оЁ 60 % сотрагеа +ю то4е15: Вапдот 4, КОММ, Г.В, АВ, апа ОТ. 

О5сизяют апа Сопсияоп5. Ц уаз юпа Фа уБеп изше АВ а|еогИБипз, отеайег ассигасу 15 асшеуеа, Баё Фе 
уабИиу оРфе ГУМ а|еогИфт 1$ ШоКег. ТВегеЮге, й сап Бе гесотлтепаеа4 аз а изеВи {001 Юг деесйпе СОУПТШ-19. 


Кеумога$: СОУТО-19 ащесйоп, с1азз1ИсаНоп, пасте 1еагите то4е|5. 


Еог сйаНоп: Вгоч Коцате Аплоз, [. Зшипоу, Мабоив Мо1зе Неппапп. Сотраг15оп оЁ тасБте [еагише то4е]$ юг 
согопауй$ рге@сНоп. Адуапсеа Епотеегие ВезеагсВ, 2022, уо1. 22, по. 1, рр. 67—75. 
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Введение. Коронавирус, также известный как СОУТО-19, впервые обнаружен в Ухане (Китай) в декабре 2019 г. 
Он представляет собой семейство вирусов, начиная от простуды и заканчивая тяжелым острым респираторным 
синдромом (ТОРС). Симптомы такого вируса схожи с симптомами простуды или сезонных заболеваний. Как и 
другие респираторные вирусы, он в основном передается воздушно-капельным путем во время кашля или 
чихания. Поэтому распознавание СОУТ-19 требует тщательного лабораторного анализа, а сокращение 
ресурсов распознавания является серьезной научной задачей. Всемирная организация здравоохранения (ВОЗ) 
11.03.2020 объявила СОУШ-19, вызванный ЗАВ5-Со\У-2, пандемией, поскольку во всем мире произошел 
экспоненциальный рост числа случаев заболеваний, а спрос на интенсивные койки и соответствующие 
структуры намного превысил существующие возможности. Первыми примерами этому являются регионы 
Италии. Бразилия зарегистрировала первый случай ЗАВ$-Со\У-2 26.02.2020. Передача вируса в этой стране 
очень быстро перешла от завезенных случаев к местным и, наконец, общинным миссиям, а федеральное 
правительство Бразилии объявило о национальной общинной передаче 20.03.2020. В штате Сан-Паулу с 
населением около 12 млн человек, где находится больница Альберта Эйнштейна, по состоянию на 23.03.2020 
зарегистрировано 477 случаев заболевания и 30 связанных с ними смертей, а 27.03.2020 имели место уже 1223 
случая СОУПШ-19 с 68 сопутствующими смертями. Для замедления распространения вируса в штате Сан-Паулу 
были введены карантин и меры социального дистанцирования. Одним из мотивов этой проблемы является тот 
факт, что в контексте обширной системы здравоохранения с возможным ограничением тестирования ЗАК5- 
СоУ\У-2 нецелесообразно тестировать каждый случай, а результаты тестов могут быть использованы при 
проверке только целевой субпопуляции. Целью работы является построение на основе машинного обучения 
модели, способной прогнозировать обнаружение ЗАК $ -СоУ-2 по медицинским данным. Для этого проводится 
сравнение различных классификационных моделей машинного обучения и определяется лучшая из них с целью 
прогнозирования коронавирусов. Сравнение основано на лицах в классе |, т.е. с положительным тестом. 
Поэтому необходимо определить модель машинного обучения с лучшим отзывом и Е1-баллом для класса 1. 
Материалы и методы. За основу принят набор данных с открытым исходным кодом из израильской больницы 
Альберта Эйнштейна в Сан-Паулу. Для исследования использованы модели машинного обучения: КапдотЕоге$$ 
(КЕ), К-ближайший сосед (КММ), Машина опорных векторов (ЗУМ), Логистическая регрессия (Г.В), Дерево 
решений (ОТ) и АдаВооз (АВ), а также 10-временная техника перекрестной проверки. Проведена оценка 
некоторых показателей производительности машинного обучения, таких как точность, отзыв и оценка Е1. 
Результаты исследования. Из 5644 человек, протестированных во время пандемии СОУТ-19, 5086 человек 
дали отрицательный результат и 558 человек — положительный. При этом поддержка машинных векторов 
показала лучшие результаты в обнаружении коронавируса с отзывом — 75% и оценкой Е1 — 60% по 
сравнению с моделями: Капдот а4гШ, КММ, ГВ, АВи ПТ. 

Обсуждение и заключение. Установлено, что при использовании алгоритмов АВ достигается большая 
точность, однако стабильность алгоритма [ЗУМ является более высокой. Поэтому его можно рекомендовать 
как полезный нструмент для выявления СОУТ-19. 


Ключевые слова: выявление СОУТО-19, классификация, модели машинного обучения. 
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1. игодисНоп 

ТВе согопауилаз 1$ а уегу зеуеге асе гезриаюгу зупдготе салзе4 Бу Фе ЗАВ$-СОУ-2 ушиз. ТЬ$ упи$, 
УТ!сВ сап шЕсё Битап$ ог апипа[5, \уаз 41зсоуеге ш фе СЫшезе гег1оп оЁ \Мирап, тоге ргес1зе!у ш Фе ргоушсе оЁ 
Нибе!, Фит» Фе рпеитота ер14епис ог Лапиагу 2020 [1,2]. [ 15 фегеюге фе зеуепй Биштап сотопауйи$. То 
еуегуопе'5 зигриве, 11$ уп$ зргеа4 моПА\лае, саязте 318,599 4еа $ апа 4,806,299 шЕеме4 регзоп$ [3]. 

ЗАК 5-СоУ-2, ЗАК $-СоУ апа МЕК$-СОУ (Мпаа Еа$ Кезриаюгу Зупаготе Согопауйл1$) сайзе зеуеге 
рпеитоша УИ а поцаШу гайе оЁ2.9 %, 9.6 % апа 36 % гезресиуейу [4—6]. 

ТЬе офег Юг уши$е$, пате!йу ОС43, МГ.63, НКУТ, апа 229Е, аге гезроп$е юг Шпеззез ге|айе4 ю шПа 
зушрюищ$ [7]. 
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[Е звош@ Бе по{е4 {вай зшсе фе Соу19-19 ер14епис, Ффеге Ваз Бееп шас зресшайоп абоиё Ше опот оЁ 1$ 
уаз [8]. Зоте за1А Фа И \уаз Ше гези оЁ ууотК 4опе ш а 1абогаюогу. Но\уеуег, айег за ез сопаис{е4 оп эепейс 
аа, 1$ Вурое$1$ \уаз 4151115364 [9]. Апа[у$15 апа сотранзоп \/ИВ е эепотез оЁ ргемоц$1у Кпо\’п согопауй$е$ 
с1еа у зВо\/ фа ЗАВ$-СОУ-2 1$ ЧИегепе Нот офег согопауйзе$ [8, 11]. Тве ушаз$ гезроп$ Ше Юг Фе согопауйа$ 
(ЗАК $-СОУ-2) 1$ зииПаг ю Ше ЗАВ$ ушиз оРБаб [2]. Твиз, е Соу19-19 ушиз 15 БеПеуеа ю Бауе опешаеа Нот а 
Ба согопаугиа$ а Бесате шесНоиз$ 0 Витапз \%ЮПе асаиите эепез зрес!Нс 10 рапгоПп согопауйа$ез. [ зВошША Бе 
пое4 Фа{ Фе асфша| сайзез оЁ Соу14-19 аге $НП ипсеаг. 

Тве зутрю$ оЁ Соу19-19 аге зппПаг ю Фозе оР зеазопа| Ни. ТВе 4еазе 15 тоге зеуеге ш фе е!4епйу апа 
ш реор!е \По аге уитегаЫе {0 сева спгошс 1зеазез. РаНеп УИ Соу14-19 сап Вауе зутрю$ гапеше Нот тИа 
То зеуеге. ТВе п105 соттоп зутрюти$ аге еуег (83 %), соизВ (82 %) апа Бгеае$зпез$ (31 %) [12]. ш райет$ \и 
рпеитошаа, Пе Х-гау оРе [ап5$ зВо\$ питегои$ поШез ап4 этоипа 21аз$ орасйу [12, 13]. 

СазфопиезИпа| зутропл$ аззостайе4 ул райеп ул Соу14-19 шс[аде уотй ше, Фагоеа, ап аб4опита1 
рат [12, 14]. 

\е а[50 ее а 4есгеазе ш [утрвосуе$ ап4 еоз1ор|ВИ$, |о\уег Ваето21о бт |еуе]з, ап4 ап шсгеазе ш \уие 
Ыоо4 сеП$ ап пеигорЬИ$ [15—18]. 

Те тапНезайноп оЁ Соу14-19 ш сЬПагеп 15 аегепе Нот Фа ш ада. ш сЬПагеп, бе зутрюпав аге шПЧ. 
Но\еуег, ш зоте сЮПагеп, \уе Вауе зееп зеуеге апа Ёайа| сазез [19—27]. 

Таке аП ошег уазез, Со\14-19 1$ Напзпие таш!Му Бу Ше гезриаюгу тоше. Атопе Фезе гоиез$ оЁ 
{гап$и11$$1оп, \е Вауе агор!ей шапзпа1$1юп, УС 15 Фе то5ё ум4езргеаа [28, 29]. ОШег тапзп$$1юп гоиёе$ ех1$, 
пате|у ше Ёеса| гоще, у1а заПуа. шдееа, ЗАВ $-СоУ-2 КМА уаз Юипа ш е $1001 оРа раме зу Соу14-19 [31]. 
ЗАК $-Со\У-2 ВМА сап Бе деесме4 оп шапита{е зигРасез (4оог Вап4]е$). Реоре \То Бауе Бееп ш сотасЕ зу Шфезе 
зиггасез$ соШа Бе сотаттаеа [29]. 

Тыз$ шо4е|! уШ таке й роззе ю 14епКу розуе ап пезайуе сазез Нот фе Чайазеё зе ап@ е 
е!етлеп{5 гезроп Ле юг СОУПШ-19. ТВе ргорозе4 рге@ сноп то4е| епзигез аё й ‘гасК$ Ше гези5 гесагИ те 115 
ерЧепис зИаайоп $0 фай фе Визе есопопис 1053е$, {фе зргеа4 оЁ Фе соттипйу, е атоцпе оЁ деасЬтеп{ зосла| 
эзеп$ сап Бе 4еесе ап4 а ргес1зе 4ес1з1юп сап а!зо Бе та4е ассог@те]у. Тв тефо@Я \уШ аПо\ воуегитег 
айфоне ю риё ш расе ргеуепнуе теазигез Базе оп очг Виаге \уогК 0 рге41с{ Фе опзе{ оЁ 15 415еазе ш фе Виеаге. 

2. Ома Везоигсез апа Ме#о4$ 

Тье Чаёазеё ‘зе \уаз ироа4е4 ю Касзе. Ш 15 ореп зошсе ап ауаЦаМе оп 5 ШК 
Касбе.сот/етует4ааи/соу1а19. ТЫ5$ Чаёазей сошашз апопупихе Чайа ш ассог4апсе \уй Безё имегпайопа| 
ргасисез ап райеп гесоттепдайоп$ ай 1е [згаеша АШек Етуешт Нозриа| ш Зо Рашо, ВтатИ. ТЬ1$ зесНоп 
Чезсгез Ше ргорозе4 арргоасй ап4 а дейаПе4 оуегуе\у оЁ Ше 1азКз. ТБезе 1азК$ сап Вер ю ипаег\апа ап@ ехгас{ 
Кпо\Ледэе Нот СОУТО 19 даа, умер сап Бер соипилез сопёашт Фе зргеаЯ оЁ е у!газ, га1зе а\агепезз, [алпсв 
шшайуез, Чеегите Ш пеайоп Ваз а роз@уе е ес ог поф, 14епй Ку офег гасюгз аЁЙесйпе Фе упиз, ес. ТЬ$ УШ 
аПо\ соипитез ю ргераге Юг \\Ва{ тау Варреп ш фе пеаг Вибаге. Т1$ сошА Вер зауе Пуез ап а|Пемайе 1е азопу. 
Ерептюое1са| шогпайоп шса4ез уамои$ свагасет$Нс$ оЁ Пе сазе зи еа, шсате сазе 14епиЙЯсайолт, асе, зех, 
{агое{ уаше, 1[утрвосу(ез, [еиКосу(ез, топосуйез, Бсо3, ес. 

2.1. аа Рге-ргосе$$ то 

ш аа апа[уз1$, Фе по5ё ппропапе $ер 1$ рге-ргосеззте. Но\еуег, И 1$ пой с1еаг \’Ваё те о4$ о рге- 
геантепе Фе апог изе4. Т№$ рагЁ плиз Бе сотр!ееа. 

2.2. Оача Тгап$Ююгтайоп 

Тве даа 15 тапзКогте4 10 Бе ргосеззе апа %юоте4 1ш. х15 ог Ашег ргосеззше. АП аза \уеге погта|те4 ю 
Бауе а теап оЁ 7его ап4 а ипй з(апдага деуайоп. \УИВ а дайазей сошатшите 111 сВагасепзс$, даа пише ейтптаеа 
1115$ше уаез (78 свагасёет1$с$) ап геате4 ппропап{ сБагас{етзНс$ (33). ТЬ1$ ехрогаюгу апа1уз1$ оЁ Пе Дафа а150 
аПо\м\ме ‘5 10 14епыРу Мо саесопез оР сБагасетзИсз, патейу упиз-ге]айе сКагас{ег15с$ ап Ыоо4-геаеа 
спагас{ет1$Ис$. ТВе {агое{ уа[ае 15 Чу14е4 шю мо сайезог!ез УмсЬ аге пезайуе сазез соде4 Бу 0 ап роз1уе сазез 
со4еа Бу 1. 

Тре а{азеё от Ше 15гаеа АбегЕ Етует НозрИа| ш Зао Рашо 1$ Чту14е4 шо фапите ап4 {е${ даа. 70 % 
О 1е аа 1$ изе4 Юг рге@сйуе то4е] ‘гашие, ап е геташте 30 % 1$ изе4 Юг 1е5Ипе. Те оБеснуе оЁ то4е| 
{таштае 15 ю а4ар! фе то4е! изте 4айа бот те гашише $е{. АЙег фе то4е! 15 Когтеа, е ргедсйоп по4е!$ зоипа 
1ез1е4 0 суашае реюогтапсе ш Ше {е51 дабазе($. 
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2.3. Тве Ргорозед Моде 5 

ТЬ$ зесйоп Чезсгез фе Чегет пасте 1еагите то4е!5 изе4 т 1$ рарег. ТБезе то4е|5 аге: Капдот 
Риз (ВЕ), К-р№аз С1озе Меебог$ (КММ), Глиеаг Зиррог Уесюг МасЬше ($УМ), Го215йс Кестез$юп (ГВ), 
Пес1$1оп Тгее (ОТ), апа АдаВооз (АВ). 

Вапдот Еоге$ (ВЕ) 

Капдот К№гез{$ (ВЕ) ог гапдот 4ес1з1оп Югез5 \меге Нг$( ргорозе4 ш 1995. ТЫ 1$ а оепега| с1азз1ИсаНоп 
фашше тефо4 фае 1еп4$ ю ууогК Бейег ап га юпа! 4ес1$1юп {ее с1азИсаНоп тейо4$ (Сапзале е# а1., 2019). 
Пес1$1оп {геез аге фе Кпдатепа ВЕ с1аз1Нег$ па{ уое Юг еасВ ог 1е Югесаз{5, ап4 1е зигутуа[ ргед1сНоп 1$ Базе4 
оп Ше та]огиу уоНпе тео ш еасВ ‘гее (Вгейталп, 2001). ТБе ассигасу оЁ еасб шее ап фе ш4ерепдепсе оЁ фе 
(еез Кот еасб офег рго\Ае Фе гепабБПцу оЁ Фе с!а551Исайоп. \\е изе4 100 теез ю ргедсё мо 1агзей с1аз$ез, 
зигу!уа| ог 4еай оРГрайепз \уИВ ВБераНй$. 

№ аге5 МеюНБог (КММ) 

Те К-Мезё МеоВБог (КММ) сИазз1ег 15 опе оЁ Ше поз соттошу зе с1аз$1ИсаНоп а1еогИбив. Т№$ 
а] еогИфит сап Бе изе4 шт зеуега] аррПсайопз. 1 зауез аП уаПА аб щез ап с1а551Йез пе\/ аби щез ассогАте 10 ет 
зииЙагиу аптепз1оп. КММ 1$ а %айзИса| гесосп оп то4е| тейфо4 юг д&есипе фе Ч егепе с1аззез оРа п1о4е]. А 
{тее Чаёа згас@ге 1$ изе4 ю деегите пе Ч15апсе Бебуееп Фе рош оЁ пиегез( ап4 Фе рош5 ш Фе цаште 4айа$ей. 
ТВе айтфще 15 сазуеа Бу Из пеебогз. ш 1е с1аз1ИЯсайоп те@о4, Фе уаше оЁ К 1$ а\ауз а розшуе ищееег 
с1озе$1 10 Ше пе1еБог. Те пеагез( у1$10п$ аге з@ес{е4 гот а $61 о# с1аззез ог ргорейу уаез ое оБ]есе. 

Зиррог Уесог Масвте (ЗУМ) 

5УМ-сопоПеа [Теагите тефо4 1$ изе4 Рог с!аз$1Йсайоп ап геотезз1юп [29]. ТЬ1$ азогибт 15 а ге!айуе1у 
пе\ арргоасй ап4 Ваз реогтеа \уеП т гесеги уеагз. ТВе ЗУМ с1азз1ег 15 Базе4 оп Ппеаг с1аз51Негз ап@ ш 1е ааа 
зерагайеа Бу а го\у, Ше ЗУМ 150[а{ез фе оБ]есё5 ш 1е зресеа с1аззез. [ сап а150 14епй Бу ап4 с!аз5у шуапсез {вай 
аге по{ зирроцеа Бу 1е даа. ТВе ошу ежепз1оп оё {1$ а!еогИфт 15 0 регогт а геотезз1оп апа[у515 ю оМат а Ппеаг 
КпсНоп, ап4 апо ег еж{епз1оп {еасНез ю с1азз1Ку ве е!етеп$ ю оМат а с1азз1Ясайоп оЁ ша!у1а1 @етеп. 

Г.0215Ис Вестезчоп Моде (Г.В) 

Го21$Яс геотеззлоп 15 Ше соггезропте геотезюоп апа[у$15 фай зВошШ4 Бе регопте4 \уВеп фе 4ереп4ете 
уама Ме 1$ Чсвоютотз (Бпагу). ГлКе аП геогез$1оп апа[узез, 10215Ис геетез$1оп 15 рге@сйуе апа[уз1$. 1 1$ изе4 0 
Чезспбе Пе дайа ап ехр!ашт Фе ге!анопзЬар Бебмееп а Череп4епе Бтагу уапае ап4 опе ог тоге поглта|, ога], 
ииегуа| ог огфта| ш4ереп4епЕ уама ез, героге [30, 31]. ТЫ$ арргоасв аззитлез аё те шагу гези ЮПо\$ а 
попа зелий оп. 

Лес15юп Тгее (ОТ) Моде 

Тве Рес1$1оп Тгее 15 а сопноПеа Теагите тео4 {ай 1$ изе4 1ю зо[уе с1аз1ИсаНоп апа геэгез$1оп ргоеитв, 
Би И 1$ тоге и5е4 1ю зо]уе с1азз1Исайоп. Тиз 15 а ро\ег  с1аз51Ясайоп те#о4 Гог Ч415еазе рге@сНоп. ТЫ$ 1$ а шее 
то4е| уВеге 1е пуегпа! поде$ гергезепЕ йе спагас{ег1$Ис$ оРа даа зе, фе БгапсВе$ гергезеп( ше 4ес1$10п пез, апа 
еасВ |еаР по4е гергезеп{5 а гезай. ТВе 4ес1з1оп тее соп$155 оЁ мо по4ез, а Честз1оп по4е ап4 а 1еаЁ по4е. Бес1з1юп 
подез Вауе тире БгапсВез ап4 аге изе4 {+0 таке а 4ес151оп, Пе |еаЁ по4ез аге 1е гезий оЁПозе 4ес1$10п$5. 

Моде! АдаВоо5( (АВ) 

АдаВоозь зВогё Юг “Адарнуе ВоозИпэ”, 15 Ше Нгзё Боо5ё а!еогИбт ргорозеа Бу Егеипа ап Зсрарие т 
1996. 5 соа| 1$ © ции \еаК рге@сюог5$ шю згопе ргед1сфогз 10 зо]уе с1аз$1ИсаНоп рго ета. Еог с1аз$ИсаНоп, ше 
Йпа] едчайоп сап Бе риё ипаег фе Беа@те Бе!о\\: 

Е(х)= мет Ут а От Г тео) (1) 
\У!Пеге [„ 4епо{е$ Ше \уеаК с1аз$1Иег т ап т 4епо{е$ Фе сотгезроп4 те \уе1е{. АдаВоо$Ё сап Бе изе4 Юг се 
тесост оп, аз И 1$ а %бапдата а|зогИфил Рог 4еесйпе №сез ш ппазез. АдаВоо3{ 1$ #23, гедитез по зер, апа 1$ зпар!е 
ап еазу №0 ргоэтат. Р1аз, и Баз пе НЙехИиу 1ю Бе аЫе ю Бе сотбте4 \иИВ апу тасНте 1еагите а]еогИбт. 

2.4. ЕуашаНоп о? Ре огтапсе Меа$иге$ 

Еог фе сотранзоп оЁ фе 41егепе с]аз;1ЙсаНоп а!еотИБлаз изе4 ш #15 рарег, зоте тей1с$ \еге еуаша(е4. 
ТВезе аге ассигасу, гесаП, апа Е1-5соге. ТВезе тей1с$ аге сайсиа{е4 Базе оп {гие розшуез (ТР), ае пегануез (ТМ), 
а5е розуез (ЕР), ап4 #а15е пегайуез (ЕМ). ТВе збап4аг41те4 сопРаз1юп тафх Шизгаез фе г@айоп$ р Бебуееп 
с1а5$1ИсаНоп гези ап ргед1се4 с1аззез. ТВе 1еуе[ оЁ {Ве с1аз$1ИсаНоп регРюгтапсе 1$ сасиае4 Бу Фе питбег оЁ 
затр[ез соггесЙу ап шсоггесЙу с1азз1Не4 ш еасВ с1азз. 


ТВе ассигасу 15 са1сшае4 Базе оп Фе {ю{а1 питбег оЁ соггесй рге@сНопз, 4еНпе4 аз ЮПо\\з: 
ТР+ТМ (2) 


Ассигасу = ————— 
% ТР+ЕМ+ТМ+ЕР 
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ВКеса|, ог зепзиуИу, 15 Ше ргорогйоп оЁ тае розуе рге@сНоп$ аё Вауе Бееп соттесЙу 14епийе4, 4ейпе@ аз 


ЮПоуз: 
Весай = —"#_— (3) 
ТРЕМ 
Тве Е1 зсоге 1$ фе Вагтотис теап оРассигасу ап4 гесаП, апа # 15 са1сщаеа Ъу: 
бсоге Е1 = Е (4) 
ТР+>(ЕР+ЕМ) 


3. Вей 

Тье об]еснуе оЁ 1$ рарег 15 0 сотраге Ше @1Шегеп по4е15 оЁ пасте 1еагише Юг Фе а&есйоп оЁ 
согопауй$. Опг базК уаз ю Япа ой \мсеН тасЬте [еаглте то4е! Ваз Фе Без гесаП апа Ё1-зсоге ог СЛазз 1. ТВе 
1еагите тасбте по4е!5 зе аге: Кадот г, К-пеагезё пелеВБог, 1091$Ис теогезз1юоп, заррой уесюг тасНше, 
АдаВооз%, ап4 4еслз1оп {гее. Ош оРа 1ю[ о# 5,644 реор!е 1е$е4 Рог СОУТО-19, 5,086 реор[е 1е5е4 пегануе апа 558 
реор1е 1е$е4 розШуе. ТВе гезаИ$ оЁ ошг за4у аге ргеземеа ш Е1еиге 1 апа Е!еаге 3. ТБезе гези$ зПо\у паЁ Пе 
уесюг-тасНше сауе Бе{ег гези{5 УЛЕЙ а гесаП оё 75 % ап4 ап Е1 $соге оЁ 60 %. Тве ЧЕгепе [еагит сигуез \еге 
а150 тасе ш ог4ег ю ип4ег%апа фе рВепотепоп оЁ оуег-Ё пе ап4 ипаег-Р то Е1оите 2. п4ее4, фе 1еагите сигуе 
15 уегу уеП Кпо\уп © 4айа 5с1епИз5, Фе 1еагите сигуе зВо\уз фе еРЯсепсу ап адиаШу оф ]еагише оЁ оиг пасте 
[еагише то4е|. Геагише сигуез аге \14еЙу изе аз а ЧаопозИс 1001] ш тасБше Пеагише Юг аеогИбил$ фай 
шсгететаПу [еагп а гашше 4аёа зе{. ТЬ1$ пеапз {фа{ \\е шсгеазе оиг Чайазе{ Бу а сецат $ер, ап4 еп \е зее Пе 
реггтапсе оЁ опг под]. ТВе то4е! сап Бе еуамае4 оп Ше таште Чабазе{ ап оп фе ехсерНоп уаПЧайоп Ча{азей 
аВег еасБ ир4же Фигте гашипе, ап4 # {гасез Бе теазиге4 регогтапсе. Т№1$ сап Бе гергезеие4 аз а сигуе. 


Капаотроге$* АдаВоо$т 
[[91 4] [91 4] 
[11 $] — [9 7]] 
ргес1$10оп  геса11 11-5соге  зирроге ргес1$10п  геса!] 11-зсоге  зиррог{ 
0 0.89 0.96 0.92 95 0 0.91 0.96 0.93 95 
1. 0.56 0.31 0.40 16 1 0.64 90.44 0.52 16 
ассигасу 0.86 111 ассигасу 0.88 111 
масго ауё 0.72 0.64 0.66 111 тмасго а\у8 0.77 0.70 0.73 111 
метей{е амё 0.84 0.86 0.85 111  ие1еей ауё 0.87 0.88 0.87 111 
км Оес151опТгее 
[88 7] [86 9] 
[8 8]] —_ [11 5]] 
ргес1$10п  геса11 +1-5соге  зиррог ргес1$10п  геса1] +1-5соге  зиррогё 
0 0.92 0.93 0.92 95 0 0.89 0.91 0.90 95 
1 0.53 0.50 0.52 16 н 0.36 0.31 0.33 16 
ассигасу 0.86 111 ассигасу 0.82 111 
тасго а\уз 0.72 0.71 0.72 111 тасго а\ув 0.62 0.61 0.61 111 
ме1ей{ей ауё 0.86 0.86 0.86 111 ме1еме4 ауё 0.81 0.82 0.81 111 
1021511с_Верге$$1оп $УМ 
[92 3] [[83 12] 
[19 6]] [4 12] 
ргес1$1оп геса11 +1-соге  зиррогф ргес1510п  геса 11-зсоге — зиррогЁ 
0 90.90 0.97 0.93 95 0 0.95 0.87 0.91 95 
1. 0.67 90.38 0.48 16 1 0.50 0.75 0.60 16 
ассигасу 9.88 111 ассигасу 0.86 111 
‚ВаеГО а\8 0.78 0.67 0.71 111 тасго ау 0.73 0.81 0.76 111 
мете[фей а\уё 0.87 0.88 0.87 111 иеуенЕей амё 0.89 0.86 0.87 111 
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Ешюлге 3 зВо\из фе регогтапсе оЁ фе @егеп( тасбше 1еагите а|зогИйил$ ассогате ю Фе регогтапсе 
теазиге$ изе4 т 15 рарег. \\е зее фа{ Гог гесаП апа Е1-зсоге, [ЗУМ ошреггиаз фе офег тасНше [еагише 1104615 
изе4, пате[у Г.В, КММ, ВЕ, АВ, апа ОТ. Еог ассигасу, ГВ 1$ тисН Бейег Пап фе оегз. Аз Юг ассигасу, \е Ипа 
Ша ГВ апа АВ регЮюгте4 Бецег ап фе офег то4е]5. ш 15 рарег, \е спозе гесаЙ апа Е1 зсоге {0 шеазиге е 
регогппапсе оЁ фе по4е1. ВесаП аПо\ме4 и {ю соттесЙу 14еп у Фе Со\19-19 розйуе {езё заб]есё$ атопз а е геа1 
роз1уе сазез. Аз Юг Ше Е1 зсоге, \е изеа # Бесаизе \е Ва ап пиба!апсе Бебуееп Ч егепе с1аззе$, 1.е., розуе апа 
пегайуе сазез. 

4. 015си5$10п апа Сопси$10п 

Тве даа изе4 ш 1$ рарег \аз соПесе4 аё {Пе 15гаеШа АЪБем Етзеш НозрИа| ш Зао Раию, Вга71. АЙег ап 
ехр!огабюгу апа[уз1$, мо саегопез оЁ сВагасет$Исз \уеге 14епийеч. ТКезе аге ше свагасегзИс$ геайе4 ю Фе уаз апа 
Фе сВагасмеп$Ис$ ге|аёе4 ю Ше Боо4. Ой оГа ю1[ оЁ 5,644 реоре 1е$е4 хип СОУТО-19, 5,086 реор[е 1е$е4 пегайуе 
ап 558 реоре {е$е4 роз@уе. Тве гези$ оЁ 515 за4у <еаЙу Шизиайе фа ш г@айоп ю ог 2оа|, пасте уесюг 
зиррог зВо\еа Бецег гезиЙ$ ш согопауйз даесНоп У\/ИВ а гесаЙ ог 75 % ап ап Е] зсоге оЁ 60 %. ТЫ$ со-са1сиайоп 
уаз 4опе уу Ше ошег пасте |еагише по4е]$, пате[у фе Кадот агШ, Ве К-пеаге$( пелеНБог, Ше 1051$ с геогез!оп, 
Ше АдаВоо$, ап Фе 4ес1з10п ее. Аз заср, 1$ тоде| сап Бе изе И Юг Фе 4а2п0$1$ оЕ СОУТ-19. Ноуеуег, И 1$ 
роз1Ые ю орйтияе Ше рагатеег$ оЁ {1$ по4е! т ог4ег ю ппргоуе Из регогтапсе. 

АЁег Фе апа[уз15 оР Фе 1еагише сигуе ш Е1оиге 2, уе Вп@ Фа арагё Нот Фе зиррогипе зепзог, офег 
тасрте еагпте то4е!$ сап Бе зе Юг фе а&есноп оЁ СОУШ-19. Твезе шса4е АдаВооз(Ё ап4 К-пеагез 
пеоНБог. шаее4, уе йп4 а! 1Ё ме регЮгт а Ше тоге а4уапсе4 орипихайоп оЁ Ве рагатеегз оЁ Фезе то4е[, Шеу 
сош Бе сапа Чаез Юг Фе Фаепоз15 оР СОУПТ-19 Бесалзе фе Ч1егепсе Бебмееп Ше Пеагише зсоге сигуе ап4 е 
уаПЧайоп зсоге сигуе \уош Вауе гедасе4 е то4е!'$ аБПиу {0 сепегаПте. 
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